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Abstract 

We describe Superlso v3.0 which is a public program for evaluation of flavor physics 
observables in the Standard Model (SM), general two-Higgs-doublet model (2HDM), 
minimal supersymmetric Standard Model (MSSM) and next to minimal supcrsymmot- 
ric Standard Model (NMSSM). Superlso v3.0, in addition to the branching ratio of 
B — ► X s 7 and the isospin asymmetry of B — > if*7, incorporates other flavor observ- 
ables such as the branching ratio of B s — > fi + pT , the branching ratio of B — » tu t , the 
branching ratio of B — ► Dtv t , the branching ratio of K — * [w^ and the branching ratios 
of D 8 — > tv t and D s — * [iv^. The calculation of the branching ratio of B — > X s j is also 
improved, as it includes NNLO Standard Model contributions. The program also com- 
putes the muon anomalous magnetic moment (g — 2). Nine sample models are included 
in the package, namely SM, 2HDM, and mSUGRA, NUHM, AMSB and GMSB for the 
MSSM, and CNMSSM, NGMSB and NNUHM for the NMSSM. Superlso uses a SUSY 
Les Houches Accord file (SLHAl or SLHA2) as input, which can be either generated 
automatically by the program via a call to external spectrum calculators, or provided 
by the user. The calculation of the observables is detailed in the Appendices, where a 
suggestion for the allowed intervals for each observable is also provided. 
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*Electronic address: mahmoudi@in2p3.fr 
URL: http://superiso.in2p3.fr 



1 



Contents 



1 Introduction 

2 Content of the Superlso v3.0 package 

2.1 Parameter structure 

2.2 Main routines 

3 Compilation and installation instructions 

4 Input and output description 

4.1 SLHA input file 

4.2 SM main program 

4.3 mSUGRA inputs 

4.4 AMSB inputs 

4.5 GMSB inputs 

4.6 NUHM inputs 

4.7 CNMSSM inputs 

4.8 NGMSB inputs 

4.9 NNUHM inputs 

4.10 2HDM inputs 

5 Results 

6 Conclusion 

Appendix A: QCD coupling 

Appendix B: Evolution of quark masses 

Appendix C: Wilson coefficients at matching scale 

C.l Standard Model contributions 

C.2 Charged Higgs contributions 

C.3 Supersymmetric contributions 



Appendix D: Renormalization group equations 

D.l RGE in the standard operator basis 

D. 2 RGE in the traditional operator basis 

Appendix E: Calculation of flavor observables 

E. l Branching ratio of B —* X s j 

E.2 Isospin asymmetry of B — ► K*j 

E.3 Branching ratio of B u —* rv T 

E.4 Branching ratio of B — > Dtv t 

E.5 Branching ratio of B s — ► 

E.5.1 Supersymmetry 

E.5. 2 2HDM 

E.6 Branching ratio of -fT — > /Ui/„ 

E. 7 Branching ratio of D s — > 

Appendix F: Muon anomalous magnetic moment 

F. l Supersymmetry 

F.2 2HDM 

Appendix G: LHA file format for 2HDM 

Appendix H: Suggested limits 

Appendix I: Useful parameters 

References 



List of Tables 



1 Yukawa couplings for the four types of 2HDM 

2 dj numbers 

3 m kli vames 

4 m ku values 

5 values 

6 Useful numbers 

7 Useful numbers for 

8 Higgs and MSSM particles mass limits .... 

9 Limits for implemented observables 

10 Input parameters 



1 Introduction 



Along with the direct searches for new physics and particles, indirect searches appear as 
very important complementary tools to explore physics beyond the Standard Model (SM). 
The presence of new particles as virtual states in processes involving only ordinary particles 
as external states provides us with the opportunities to study the indirect effects of new 
physics. This is the case for example of rare B decays or the anomalous magnetic moment 
of the charged leptons. We investigate here the extension of the Higgs sector in the two- 
Higgs-doublet model (2HDM) and supersymmetry (SUSY), which are the most popular 
extensions of the Standard Model. 

Phenomenological interests of indirect constraints are numerous. This includes, in first 
place, tests of the SM predictions. By studying the indirect constraints it is also possible to 
reveal indirect effects of new physics. This aspect can be complementary to direct searches 
or even used as guideline. Finally, the results of indirect searches can be employed in order 
to check consistencies with the direct search results. 

The main purpose of the Superlso program is to offer the possibility to evaluate the most 
important indirect observables and constraints. 

Superlso [1] was in its first versions devoted to the calculation of the isospin symme- 
try breaking in B — > K*j decays in the minimal supersymmetric extension of the Standard 
Model (MSSM) with minimal flavor violation. This observable imposes stringent constraints 
on supersymmetric models [2], which justifies a dedicated program. The calculation of the 
b — > S7 branching ratio was also included in the first version and has been improved by 
adding the Standard Model NNLO contributions since version 2. Also, a broader set of 
flavor physics observables has been implemented. This includes the branching ratio of 
B s — ► [i + fi~, the branching ratio of B u — ► tv t , the branching ratios of B — ► D°tu t and 
B — ► D°eu e , the branching ratio of K — > \iv„ and the branching ratios of D s — > tv t and 
D s — * [iv The calculation of the anomalous magnetic moment of the muon is also imple- 
mented in the program. 

Superlso has been extended to the general two-Higgs-doublet model since version 2.6 and 
to the next to minimal supersymmetric extension of the Standard Model (NMSSM) since 
version 3.0 [3]. 

Superlso uses a SUSY Les Houches Accord file [4,5] as input, which can be either generated 
automatically by the program via a call to SOFTSUSY [6], ISAJET [7], NMSSMTools [8], or 
provided by the user. Superlso can also use the LHA inspired format for the 2HDM, gener- 
ated by 2HDMC [9], which will be described in Appendix[G} The program is able to perform 
the calculations automatically for different types of 2HDM (I— IV), for different supersym- 
metry breaking scenarios, such as the minimal Supergravity Grand Unification (mSUGRA) 
also called Constrained MSSM (CMSSM), the Non-Universal Higgs Mass model (NUHM), 
the Anomaly Mediated Supersymmetry Breaking scenario (AMSB) and the Gauge Medi- 
ated Supersymmetry Breaking scenario (GMSB), and for the NMSSM scenarios namely 
CNMSSM, NGMSB and NNUHM. 
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In the following, we first discuss the content of the Super I so package and give the list 
of the main routines. Then we describe the procedure to use Superlso, and introduce the 
inputs and outputs of the program. Finally, we present some examples of results obtained 
with Superlso. In the Appendices, a complete description of the formulas used to calculate 
the different observables in Superlso is given for reference, as well as suggestions for the 
allowed intervals. 



2 Content of the Superlso v3.0 package 

Superlso is a C program respecting the C99 standard, devoted to the calculation of the 
most constraining flavor physics observables. Ten main programs are provided in the pack- 
age, but the users are also invited to write their own main programs, slha.c can scan 
files written following the SUSY Les Houches Accord formats, and calculates the corre- 
sponding observables. sm.c provides the values of the observables in the Standard Model 
while thdm . c computes the observables in 2HDM types I-IV and requires 2HDMC [9] for 
the generation of the input file containing the Higgs masses and couplings. The main pro- 
grams msugra. c, amsb. c, gmsb. c, and nuhm. c have to be linked to the SOFTSUSY [6] and/or 
the ISASUGRA/ISAJET [7] packages, in order to compute supersymmetric mass spectra and 
couplings within respectively the mSUGRA, AMSB, GMSB or NUHM scenarios. The pro- 
grams cnmssm.c, ngmsb.c, and nnuhm.c have to be linked to the NMSSMTools [8] program 
to compute supersymmetric mass spectra and couplings within respectively the CNMSSM, 
NGMSB or NNUHM scenarios. 

The computation of the different observables in Superlso proceeds following three main 
steps: 

• Generation of a SLHA file with ISAJET, SOFTSUSY, NMSSMTools or 2HDMC (or supply 
of a SLHA file by the user) , 

• Scan of the SLHA file, 

• Calculation of the observables. 

The last point hides a complex procedure: to compute the inclusive branching ratio of 
b — ► 57 for example, Superlso needs first to compute the Wilson coefficients at matching 
scale, and then to evolve them using the Renormalization Group Equations (RGE) to a 
lower scale, before using them to compute the branching ratio. 

2.1 Parameter structure 

The package Superlso v3 . has been deeply revised and expanded as compared to the first 
release, and many of the original routines have been modified. The code still relies on the 
definition of a structure in src/ include . h, which has been also expanded, in particular in 
order to support the new parameters of the SLHA2 format. This structure is defined as 
follows: 



6 



typedef struct parameters 

/* structure containing all the scanned parameters from the SLHA file */ 
{ 

int SM; 

int model; /* mSUGRA = 1, GMSB = 2, AMSB = 3 */ 
int generator; /* ISAJET = 1, SOFTSUSY = 2 */ 

double Q; /* Qmax ; default = M_EWSB = sqrt (m_stopl*mstop2) */ 

double m0,ml2,tan_beta,sign_mu,A0; /* mSUGRA parameters */ 

double Lambda, Mmess,N5,cgrav,m32; /* AMSB, GMSB parameters */ 

double mass_Z,mass_W,mass_b,mass_top_pole,mass_tau_pole; /* SM parameters */ 

double inv_alpha_em, alphas_MZ , alpha, Gfermi ,GAUGE_Q ; /* SM parameters */ 

double charg_Umix [3] [3] , charg_Vmix [3] [3] ,stop_mix[3] [3] ,sbot_mix[3] [3] , 

stau_mix[3] [3] ,neut_mix[6] [6] ,mass_neut [6] ; /* mass mixing matrices */ 

double Min , Ml _Min , M2_Min , M3_Min , At_Min , Ab_Min , Atau_Min , M2H1 _Min , M2H2_Min , 

mu_Min,M2A_Min,tb_Min,mA_Min; /* optional input parameters at scale Min */ 

double MeL_Min , MmuL_Min , MtauL_Min , MeR_Min , MmuR_Min , 

MtauR_Min; /* optional input parameters at scale Min */ 

double MqLl_Min , MqL2_Min , MqL3_Min , MuR_Min , McR_Min , MtR_Min , 

MdR_Min,MsR_Min,MbR_Min; /* optional input parameters at scale Min */ 

double N51 ,N52,N53,M2H1_Q,M2H2_Q; /* optional input parameters (N51...3: GMSB) */ 

double mass_d ,mass_u,mass_s ,mass_c ,mass_t ,mass_e ,mass_nue ,mass_mu, 

mass_num,mass_tau,mass_nut ; /* SM masses */ 

double mass_gluon,mass_photon,mass_ZO; /* SM masses */ 

double mass_hO ,mass_H0 ,mass_A0 ,mass_H ,mass_dnl ,mass_upl ,mass_stl ,mass_chl , 
mass_bl ,mass_tl ; /* Higgs & superparticle masses */ 

double mass_el ,mass_nuel ,mass_mul ,mass_numl ,mass_taul ,mass_nutl ,mass_gluino , 
mass_chal ,mass_cha2 ; /* superparticle masses */ 

double mass_dnr ,mass_upr ,mass_str ,mass_chr ,mass_b2 ,mass_t2 ,mass_er ,mass_mur , 

mass_tau2; /* superparticle masses */ 

doubl e mas s _nuer , mas s _numr , mas s _nut r , mas s _gr avit on , 

mass_gravitino ; /* superparticle masses */ 

double gp,g2,g3,YU_Q,yut [4] ,YD_Q,yub[4] ,YE_Q,yutau[4] ; /* gauge & Yukawa couplings */ 
double HMIX_Q ,mu_Q , tanb_GUT , Higgs_VEV ,mA2_Q , MSOFT.Q , M1_Q ,M2_Q , 
M3_Q; /* parameters at scale Q */ 

double MeL_Q , MmuL_Q , MtauL_Q , MeR_Q , MmuR_Q , MtauR_Q , MqLl_Q , MqL2_Q , MqL3_Q , MuR_Q , McR_Q , 
MtR_Q , MdR_Q , MsR_Q , MbR_Q ; /* masses at scale Q */ 

double AU_Q,A_u,A_c,A_t,AD_Q,A_d,A_s,A_b,AE_Q,A_e,A_mu,A_tau; /*trilinear couplings*/ 
/* SLHA2 */ 

int NMSSM , Rparity , CPviolat ion, Flavor ; 

double mass_nutau2 ,mass_e2 ,mass_nue2 ,mass_mu2 ,mass_numu2 ,mass_d2 ,mass_u2 ,mass_s2 , 
mass_c2; 

double CKM_lambda, CKM_A , CKM_rho , CKM_eta ; 

double PMNS_thetal2 , PMNS_theta23 , PMNS_thetal3 , PMNS_deltal3 , PMNS_alphal , PMNS_alpha2 ; 
double lambdaNMSSM_Min , kappaNMSSM_Min , AlambdaNMSSM_Min , AkappaNMSSM_Min , 
lambdaSNMSSM_Min , xiFNMSSM_Min , xiSNMSSM_Min , mupNMSSM_Min , mSp2NMSSM_Min , mS2NMSSM_Min , 
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mass_H03 , mass_A02 , NMSSMRUN_Q , lambdaNMSSM , kappaNMSSM , AlambdaNMSSM , AkappaNMSSM , 
lambdaSNMSSM ) xiFNMSSM,xiSNMSSM,mupNMSSM,mSp2NMSSM,mS2NMSSM; /* NMSSM parameters */ 
double PMNSU_Q , CKM_Q , MSE2_Q , MSU2_Q , MSD2_Q , MSL2_Q , MSQ2_Q , TU_Q , TD_Q , TE_Q ; 
int NMSSMcoll,NMSSMtheory,NMSSMupslS,NMSSMetablS; /* NMSSM exclusions */ 

double CKM[4] [4] ; /* CKM matrix */ 

double H0_mix[4] [4] ,A0_mix[4] [4] ; /* Higgs mixing matrices */ 

double sU_mix[7] [7] ,sD_mix[7] [7] , sE_mix[7] [7] , sNU_mix[4] [4] ; /* mixing matrices */ 

double sCKM_msq2 [4] [4] , sCKM_msl2 [4] [4] , sCKM_msd2 [4] [4] , sCKM_msu2 [4] [4] , 

sCKM_mse2 [4] [4] ; /* super CKM matrices */ 

double PMNS_U[4] [4] ; /* PMNS mixing matrices */ 

double TU[4] [4] ,TD[4] [4] ,TE[4] [4] ; /* trilinear couplings */ 

/* non-SLHA */ 

double mass_b_lS,mass_b_pole ,mtmt ; 
double Lambda5; /* Lambda QCD */ 

/* Flavor constants */ 

double f _B,f _Bs,f _Ds,m_B,m_Bs,m_Ds,m_K,m_Kstar,m_D,life_B,lif e_Bs , lif e_Ds ; 
/* Higgs widths */ 

double width_hO , width_H0 , width_AO , width_H , width_H03 , width_A02 ; 

/* 2HDM */ 

int THDM_model ; 

double lambda_u[4] [4] ,lambda_d[4] [4] ,lambda_l[4] [4] ; 
} 

parameters ; 

This structure contains all the important parameters and is called by most of the main 
functions in the program. 

2.2 Main routines 

The main routines of Super I so are explained in the following. The complete list of imple- 
mented routines can be found in src/include .h. 

• void Init_param(struct parameters* param) 

This function initializes the param structure, setting all the parameters to 0, apart 
from the SM masses and couplings, which receive the values given in the PDG [10]. 

• int Les_Houches_Reader (char name [] , struct parameters* param) 

This routine reads the SLHA file whose name is contained in name, and put all the read 
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parameters in the structure param. This function has been updated to the SLHA2 
format. This routine can also read the LHA inspired format for the 2HDM described 
in Appendix [Gj A negative value for param->model indicates a problem in reading 
the SLHA file, or a model not yet included in Superlso (such as -R-parity breaking 
models). In this case, Les_Houches_Reader returns 0, otherwise 1. 

• int test_slha(char name [] ) 

This routine checks if the SLHA file whose name is contained in name is valid, and if 
so returns 1. If not, -1 means that in the SLHA generator the computation did not 
succeed (e.g. because of tachyonic particles), -2 means that the considered model is 
not currently implemented in Superlso, and -3 that the provided file is either not in 
the SLHA format, or some important elements are missing. 

• int sof tsusy_sugra(double mO, double ml2, double tanb, double AO, 
double sgnmu, double mtop, double mbot, double alphas_mz, char name [] ) 

• int isajet_sugra(double mO, double ml2, double tanb, double AO, 
double sgnmu, double mtop, char name [] ) 

• int sof tsusy_gmsb (double Lambda, double Mmess, double tanb, int N5, 

double cGrav, double sgnmu, double mtop, double mbot, double alphas_mz, 
char name [] ) 

• int sof tsusy_amsb (double mO, double m32, double tanb, double sgnmu, 
double mtop, double mbot, double alphas_mz, char name [] ) 

• int sof tsusy_nuhm(double mO, double ml2, double tanb, double AO, double mu, 
double mA, double mtop, double mbot, double alphas_mz, char name [] ) 

The above routines call SOFTSUSY or ISAJET to compute the mass spectrum cor- 
responding to the input parameters (more details are given in the next sections), and 
return a SLHA file whose name has to be specified in the string name. 

• int thdmc_types (double 11, double 12, double 13, double 14, double 15, 
double 16, double 17, double ml2_2, double tanb, int type, char name [] ) 

This routine calls 2HDMC to compute the masses and couplings corresponding to the 
2HDM input parameters, and returns a LHA inspired file whose name has to be spec- 
ified in the string name. 

• int nmssmtools_cnmssm(double mO, double ml2, double tanb, double AO, 
double lambda, double AK, double sgnmu, double mtop, double mbot, 
double alphas_mz, char name [] ) 
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• int nmssmtools_nnuhni(double mO, double ml2, double tanb, double AO, 
double MHDGUT, double MHUGUT, double lambda, double AK, double sgnmu, 
double mtop, double mbot, double alphas_mz, char name [] ) 

• int nmssmtools_ngmsb (double Lambda, double Mmess, double tanb, int N5, 
double lambda, double AK, double Del_h, double sgnmu, double mtop, 
double mbot, double alphas_mz, char name [] ) 

The above routines call NMSSMTools to compute the mass spectrum corresponding 
to the input parameters (more details are given in the next sections) and return the 
SLHA2 file name. 

• double alphas_running (double Q, double mtop, double mbot, 
struct parameters* param) 

This function computes the strong coupling constant at the energy scale Q using the 
parameters in param, provided the top quark mass mtop and bottom quark mass mbot 
used for the matching between the scales corresponding to different flavor numbers are 
specified. The main formula for calculating a s is given in Appendix[A] The prototype 
of this function has changed since version 2. 

• double running_mass (double quark_mass, double Qinit, double Qfin, double mtop, 
double mbot, struct parameters* param) 

This function calculates the running quark mass at the scale Qfin, for a quark of 
mass quark_mass at the scale Qinit using the structure param, knowing the match- 
ing scales mtop and mbot. A description of the quark mass calculations can be found 
in Appendix iBl 

• void CW_calculator (double C0w[] , double Clw[] , double C2w[] , double mu_W, 
struct parameters* param) 

• void C_calculator_basel (double C0w[] , double Clw[] , double C2w[] , double mu_W, 
double C0b[], double Clb[], double C2b[], double mu, struct parameters* param) 

• void C_calculator_base2 (double C0w[] , double Clw[] , double mu_W, double C0b[], 
double Clb[], double mu, struct parameters* param) 

These three routines replace C_calculator in the first version. 

The procedure CW_calculator computes the LO contributions to the Wilson coeffi- 
cients C0w[], the NLO contributions Clw[] and the NNLO contributions C2w[] at 
the matching scale mu_W, using the parameters of param, as described in Appendix [Cl 

C_calculator_basel evolves the LO, NLO and NNLO Wilson coefficients C0w[], 
Clw[] , C2w[] initially at scale mu_W towards C0b[] , Clb[] , C2b[] at scale mu, in the 
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standard operator basis described in Appendix lD.il 

C_calculator_base2 evolves the LO and NLO Wilson coefficients C0w[] , Clw[] ini- 
tially at scale mu_W towards COb [] , Clb [] at scale mu, in the traditional operator basis 
described in Appendix ID. 21 

• double bsgamma (double CO [] , double CI [] , double C2 [] , double mu, double mu_W, 
struct parameters* param) 

This function has replaced the previous calculation of b — ► S7, which was performed 
at NLO accuracy. Here, knowing the LO, NLO and NNLO Wilson coefficients CO [] , 
CI [] , C2 [] at scale mu, and given the matching scale mu_W this procedure computes 
the inclusive branching ratio of b — > 57 at NNLO, as described in Appendix IE.1L 
The container routine double bsgamma_calculator (char name[]), in which name 
contains the name of a SLHA file, automatizes the whole calculation, as it first calls 
Init_param and Les_Houches_Reader, then CW_calculator and C_calculator_basel, 
and finally bsgamma. 

• double deltaO (double CO [], double C0_spec [] .double CI [], double Cl_spec[], 
struct parameters* param, double mub, double muspec, double lambda_h) 

This function computes the isospin asymmetry as described in Appendix IE. 21 us- 
ing both the LO and NLO parts of the Wilson coefficients at scale mub (CO [] and 
CI []), and at the spectator scale muspec (C0_spec[] and Cl_spec[]), with the ad- 
ditional input A^ in GeV. Compared to the first version, the calculation has been 
improved, and all the involved integrals have been coded in separate routines. Again, 
an automatic container routine which only needs the name of a SLHA file is provided: 
double deltaO_calculator (char name[]). 

• int excluded_mass_calculator (char name [] ) 

This routine, with the name of the SLHA file in the argument, checks whether the 
parameter space point is excluded by the collider constraints on the particle masses 
and if so returns 1 . The implemented mass limits are given in Appendix [Hi and can 
be updated by the users in src/excluded_masses . c . These limits are valid only in 
the MSSM. 

• int NMSSM_collider_excluded(char name [] ) 

• int NMSSM_theory_excluded(char name [] ) 

These two routines only apply to the SLHA file name generated by NMSSMTools, 
as they need NMSSMTools specific outputs. They respectively check if a parameter 
space point is excluded by collider constraints [11] or by theoretical constraints (such 
as unphysical global minimum). The output 1 means that the point is excluded. 

• int charged_LSP_calculator (char name [] ) 

This routine, with the name of the SLHA file in the argument, checks whether the 
LSP is charged, and if so returns 1. 
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The following routines encode newly implemented observables: 

• double Btaunu_calculator (char name [] ) 
double RBtaunu_calculator (char name [] ) 

These two routines compute respectively the branching ratio of the leptonic decay 
B u — * rv T and the ratio BR(i? u — ► rz^ T )/BR(i? u — ► ti> t ) sm as described in Ap- 
pendix IE.3I These leptonic decays occur at tree level, and we consider also higher 
order SUSY corrections to the Yukawa coupling. These routines only need the name 
of a SLHA file as input. 

• double BDtaunu_calculator (char name [] ) 
double BDtaunu_BDenu_calculator (char name [] ) 

These routines compute respectively the branching ratio of the semileptonic decay 
B -> D q tv t and the ratio BR(B -> D°rz; r )/BR(B D°e^ c ) as described in Ap- 
pendix lE~4l These semileptonic decays occur at tree level, and we consider also higher 
order SUSY corrections to the Yukawa coupling. 

• double Bsmumu_calculator (char name [] ) 

This function computes the branching ratio of the rare decay B s — * fi + fi~ consid- 
ering full one loop corrections as described in Appendix IE. 51 

• double Kmunu_pimunu_calculator (char name [] ) 
double R123_calculator (char name [] ) 

These functions compute respectively the ratio BR(K — ► ^f M )/BR(7r — * [iv^ and the 
observable R123 as described in Appendix IE, 61 These leptonic decays occur at tree 
level, and we consider also higher order SUSY corrections to the Yukawa coupling. 

• double Dstaunu_calculator (char name [] ) 
double Dsmunu_calculator (char name [] ) 

These routines compute respectively the branching ratios of the leptonic decays D s — * 
tv t and D s — ► \iv^ as described in Appendix IE. 71 These leptonic decays occur at tree 
level, and we consider also higher order SUSY corrections to the Yukawa coupling. 

• double muon_gm2_calculator (char name [] ) 

This routine computes the muon anomalous magnetic moment (5a^), considering 
one loop contributions as well as the main two loop contributions, as described in 
Appendix [FJ 
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3 Compilation and installation instructions 



The main structure of the Superlso package is unchanged since the first version, and the 
spirit of the program relies on the idea of simplicity of use. 

The Superlso package 1 can be downloaded from: 



http : / / superiso . in2p3.fr 



The following main directory is created after unpacking: 
superiso_vX . X 

It contains the src/ directory, in which all the source files can be found. The main directory 
contains also a Makefile, a README, ten sample main programs (sm.c, thdm.c, msugra.c, 
amsb.c, gmsb.c, nuhm.c, cnmssm.c, ngmsb.c, nnuhm.c and slha.c) and one example of 
input file in the SUSY Les Houches Accord format (example . lha). The compilation op- 
tions should be defined in the Makefile, as well as the path to the ISAJET isasugra.x 
and SOFTSUSY softpoint.x executable files, the path to NMSSMTools main directory, 
and the path to 2HDMC directory, when needed. 

Superlso is written for a C compiler respecting the C99 standard. In particular, it has been 
tested successfully with the GNU C Compiler and the Intel C Compiler on Linux and Mac 
32-bits or 64-bits machines, and with the latest versions of SOFTSUSY, ISAJET, NMSSMTools 
and 2HDMC. Additional information can be found in the README file. 

To compile the library, type 

make 

This creates the libisospin. a in src/. Then, to compile one of the ten programs provided 
in the main directory, type 

make name or make name . c 

where name can be sm, thdm, msugra, amsb, gmsb, nuhm, cnmssm, ngmsb, nnuhm or slha. 
This generates an executable program with the .x extension. Note that sm and slha do 
not need any additional program, but msugra, amsb, gmsb and nuhm need either ISAJET or 
SOFTSUSY, cnmssm, ngmsb and nnuhm require NMSSMTools and 2HDMC is necessary for thdm. 

sm.x calculates the different observables described in the Appendices in the Standard Model, 
using the parameters of Table [lOj 

slha.x calculates the different observables using the parameters contained in the SLHA 
file whose name has to be passed as input parameter. 

An alternative package including the calculation of the relic density, Superlso Relic [12], is also avail- 
able at: http://superiso.in2p3.fr/relic . 
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thdm.x calculates the observables in 2HDM (general model or types I-IV), starting first 
by calculating the mass spectrum and couplings thanks to 2HDMC. 

amsb.x, gmsb.x, msugra.x and nuhm.x compute the observables, first by calculating the 
mass spectrum and couplings thanks to ISA JET (for msugra. x only) and/or SOFTSUSY within 
respectively the AMSB, GMSB, mSUGRA or NUHM parameter spaces. 

cnmssm.x, ngmsb.x or nnuhm.x compute the observables, after calculating the mass spec- 
trum and couplings thanks to NMSSMTools within respectively the CNMSSM, NGMSB or 
NNUHM parameter spaces. 

For all these programs, arguments referring to the usual input parameters have to be passed 
to the program. If not, a message will describe which parameters have to be specified. 

4 Input and output description 

The inputs of slha, msugra, amsb, gmsb and nuhm are the same as in the previous version, 
but the outputs are more numerous. 

4.1 SLHA input file 

The program slha.x reads the needed parameters in the given input SLHA file and calcu- 
lates the observables. For example, the command 

./slha.x example. lha 

returns 

delta0=8.259e-02 
BR_bsgamma=2 . 974e-04 
BR_Btaunu=l . 097e-04 
Rtaunu=9.966e-01 
BR_Kmunu/BR_pimunu=6 . 454e-01 
R123=1.000e+00 
BR_BDtaunu=6 . 966e-03 
BR_BDtaunu/BR_BDenu=2 . 973e-01 
BR_Bsmumu=3 . 470e-09 
BR_Dstaunu=4 . 813e-02 
BR_Dsmunu=4 . 969e-03 
a_muon=2.045e-10 
excluded_mass=0 

where deltaO refers to the isospin symmetry breaking in B — ► K*j decays, BR_bsgamma 
the branching ratio of B — ► X s j, BR_Btaunu the branching ratio of B u — > tv t , Rtaunu the 
normalized ratio to the SM value, BR_Kmunu/BR_pimunu the ratio BR(K — > \w^) /BR(7r — > 
fii/p), R123 the ratio i?^23> BR_BDtaunu the branching ratio of B — > D°tv t , BR_BDtaunu/BR_BDenu 
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the ratio BR(i? — > Z?°r^ r )/BR(i? — > D°ev e ), BR_Bsmumu the branching ratio of B s — > fi + fi~ , 
BR_Dstaunu and BR_Dsmunu the branching ratios of Z) s — ► rz/ T and D s — ► //i/^ respectively, 
a_muon the deviation in the anomalous magnetic moment of the muon, and excluded_mass 
indicates that the point is already excluded by the direct searches if it is equal to 1. More 
details on the definitions and calculations of these observables are given in the appendices. 
If the SLHA file provided to slha.x is inconsistent, a message will be displayed: 

• Invalid point means that the SLHA generator had not succeeded in generating the 
mass spectrum (e.g. due to the presence of tachyonic particles). 

• Model not yet implemented means that the SLHA file is intended for a model not 
implemented in Superlso, such as i?-parity violating models. 

• Invalid SLHA file means that the SLHA file is invalid and misses important pa- 
rameters. 

4.2 SM main program 

The program sm . x is a standalone program which computes the different observables in the 
Standard Model. No argument is necessary for this program. The command 

. /sm . x 

returns 

delta0=8.203e-02 
BR_bsgamma=3 . 061e-04 
BR_Btaunu=l . 102e-04 
Rtaunu=l . 000e+00 
BR_Kmunu/BR_pimunu=6 . 454e-01 
R123=1.000e+00 
BR_BDtaunu=6 . 980e-03 
BR_BDtaunu/BR_BDenu=2 . 975e-01 
BR_Bsmumu=3 . 395e-09 
BR_Dstaunu=4 . 822e-02 
BR_Dsmunu=4 . 975e-03 

4.3 mSUGRA inputs 

The program msugra . x computes the observables in the mSUGRA parameter space, using 
ISA JET and/or SOFTSUSY to generate the mass spectra. If only one of these generators is 
available the corresponding #define in msugra. c has to be commented. The necessary 
arguments to this program are: 

• m,Q: universal scalar mass at GUT scale, 

• mu2- universal gaugino mass at GUT scale, 

• Aq: trilinear soft breaking parameter at GUT scale, 

• tan/3: ratio of the two Higgs vacuum expectation values. 
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Optional arguments can also be given: 



• sign(n): sign of Higgsino mass term, positive by default, 

• mf° le : top quark pole mass, by default 172.4 GeV, 

• fnb(rnb): scale independent b-quark mass, by default 4.2 GeV (option only available 
for SOFTSUSY), 

• a s (Mz): strong coupling constant at the Z-boson mass, by default 0.1176 (option 
only available for SOFTSUSY). 

If the arguments are not specified, a message will describe the needed parameters in the 
correct order. 

With SOFTSUSY 3.0.11 and ISAJET 7.79, running the program with: 
./msugra.x 500 500 50 

returns 

SLHA file generated by SOFTSUSY 
delta0=9.967e-02 
BR_bsgamma=2 . 223e-04 
BR_Btaunu=6 . 507e-05 
Rtaunu=5.905e-01 
BR_Kmunu/BR_pimunu=6 . 429e-01 
R123=9.980e-01 
BR_BDtaunu=6 . 289e-03 
BR_BDtaunu/BR_BDenu=2 . 681e-01 
BR_Bsmumu=4 . 007e-08 
BR_Dstaunu=4 . 794e-02 
BR_Dsmunu=4 . 951e-03 
a_muon=2 . 049e-09 
charged_LSP=0 
excluded_mass=0 

SLHA file generated by ISAJET 
delta0=1.009e-01 
BR_bsgamma=2 . 169e-04 
BR_Btaunu=7 . 245e-05 
Rtaunu=6.582e-01 
BR_Kmunu/BR_pimunu=6 . 434e-01 
R123=9.984e-01 
BR_BDtaunu=6 . 404e-03 
BR_BDtaunu/BR_BDenu=2 . 733e-01 
BR_Bsmumu=3 . 335e-08 
BR_Dstaunu=4 . 793e-02 
BR_Dsmunu=4 . 950e-03 
a_muon=2 . 084e-09 
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charged_LSP=0 
excluded_mass=0 

corresponding to the observables described in the previous section. 

4.4 AMSB inputs 

The program amsb . x computes the observables using the corresponding parameters gener- 
ated by SOFTSUSY in the AMSB scenario. The necessary arguments to this program are: 

• tuq: universal scalar mass at GUT scale, 

• "1-3/2 : gravitino mass at GUT scale, 

• tan/3: ratio of the two Higgs vacuum expectation values. 

Optional arguments are the same as for mSUGRA. If the input parameters are absent, a 
message will ask for them. 

With SOFTSUSY 3.0.11, running the program with: 
./amsb.x 500 5000 5 -1 

returns 

delta0=7.671e-02 
BR_bsgamma=3 . 295e-04 
BR_Btaunu=l . 097e-04 
Rtaunu=9.953e-01 
BR_Kmunu/BR_pimunu=6 . 454e-01 
R123=1.000e+00 
BR_BDtaunu=6 . 971e-03 
BR_BDtaunu/BR_BDenu=2 . 972e-01 
BR_Bsmumu=3 . 349e-09 
BR_Dstaunu=4 . 818e-02 
BR_Dsmunu=4 . 975e-03 
a_muon=-6.598e-10 
excluded_mass=l 

4.5 GMSB inputs 

The program gmsb . x computes the observables using the GMSB parameters generated by 
SOFTSUSY. The necessary arguments to this program are: 

• A: scale of the SUSY breaking in GeV (usually 10000-100000 GeV), 

• M mess : messenger mass scale (> A), 

• A%: equivalent number of 5 + 5 messenger fields, 

• tan/3: ratio of the two Higgs vacuum expectation values. 

Optional arguments are the same as for mSUGRA, with an additional one: 
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c-Grav (> 1) : ratio of the gravitino mass to its value for a breaking scale A, 1 by 
default. 



Again, in the case of lack of arguments, a message will be displayed. 

With SOFTSUSY 3.0.11, running the program with: 
./gmsb.x 2e4 5e6 1 10 

returns 

delta0=6.339e-02 
BR_bsgamma=4 . 260e-04 
BR_Btaunu=8 . 077e-05 
Rtaunu=7.330e-01 
BR_Kmunu/BR_pimunu=6 . 439e-01 
R123=9.988e-01 
BR_BDtaunu=6 . 542e-03 
BR_BDtaunu/BR_BDenu=2 . 789e-01 
BR_Bsmumu=3 . 805e-09 
BR_Dstaunu=4 . 806e-02 
BR_Dsmunu=4 . 962e-03 
a_muon=3 . 774e-08 
excluded_mass=l 



4.6 NUHM inputs 

The program nuhm . x computes the observables using the NUHM parameters generated by 
SOFTSUSY. The necessary arguments to this program are the same as for mSUGRA, with 
two additional ones, the values of [i and tua- 

• mo: universal scalar mass at GUT scale, 

• m 1 / 2 - universal gaugino mass at GUT scale, 

• A$: trilinear soft breaking parameter at GUT scale, 

• tan (3: ratio of the two Higgs vacuum expectation values, 

• //:// parameter, 

• rriA'- CP-odd Higgs mass. 

Optional arguments can also be given: 

• m v t ole : top quark pole mass, by default 172.4 GeV, 

• m^{m^)-. scale independent b-quark mass, by default 4.2 GeV, 

• a s (Mz): strong coupling constant at the Z-boson mass, by default 0.1176. 

In the absence of arguments, a message will be shown. 
With SOFTSUSY 3.0.11, running the program with: 
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./nuhm.x 500 500 50 500 500 



returns 

delta0=1.066e-01 
BR_bsgamma=2 . 014e-04 
BR_Btaunu=6 . 576e-05 
Rtaunu=5.967e-01 
BR_Kmunu/BR_pimunu=6 . 429e-01 
R123=9.981e-01 
BR_BDtaunu=6 . 300e-03 
BR_BDtaunu/BR_BDenu=2 . 686e-01 
BR_Bsmumu=3 . 650e-08 
BR_Dstaunu=4 . 795e-02 
BR_Dsmunu=4 . 951e-03 
a_muon=2.244e-09 
excluded_mass=0 

4.7 CNMSSM inputs 

The program cnmssm. x computes the observables using the CNMSSM parameters generated 
by NMSSMTools 2 . The necessary arguments to this program are: 

• mo: universal scalar mass at GUT scale, 

• m-i/2'- universal gaugino mass at GUT scale, 

• A$: trilinear soft breaking parameter at GUT scale, 

• tan/3: ratio of the two Higgs vacuum expectation values, 

• A: cubic Higgs coupling. 

Optional arguments can also be given: 

• sign(n): sign of Higgsino mass term, positive by default, 

• A K : trilinear soft breaking parameter at GUT scale, by default A K = Aq, 

• mf° le : top quark pole mass, by default 172.4 GeV, 

• rfTbirnb): scale independent b-quark mass, by default 4.2 GeV, 

• a s (Mz)- strong coupling constant at the Z-boson mass, by default 0.1176. 

In the absence of arguments, a message will be shown. 

With NMSSMTools 2.2.0, running the program with: 
. /cnmssm. x 500 500 50 0.01 

returns 

2 As the soft singlet mass m| and the singlet self coupling k are both determined in terms of the other 
parameters through the minimization equations of the Higgs potential in NMSSMTools, what we call CNMSSM 
here is a partially constrained NMSSM. 
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SLHA file generated by NMSSMtools 

delta0=9.945e-02 

BR_bsgamma=2 . 232e-04 

BR_Btaunu=6 . 478e-05 

Rtaunu=5.878e-01 

BR_Kmunu/BR_pimunu=6 . 429e-01 

R123=9.980e-01 

BR_BDtaunu=6 . 284e-03 

BR_BDtaunu/BR_BDenu=2 . 679e-01 

BR_Bsmumu=3 . 909e-08 

BR_Dstaunu=4 . 794e-02 

BR_Dsmunu=4 . 951e-03 

a_muon=2. 120e-09 

charged_LSP=0 

collider_excluded=0 

theory_excluded=0 

4.8 NGMSB inputs 

The program ngmsb . x computes the observables using the NGMSB parameters generated 
by NMSSMTools. The necessary arguments to this program are: 

• A: scale of the SUSY breaking in GeV (usually 10000-100000 GeV), 

• M mess : messenger mass scale (> A), 

• A5: equivalent number of 5 + 5 messenger fields, 

• tan/3: ratio of the two Higgs vacuum expectation values, 

• A: cubic Higgs coupling. 

Optional arguments are the same as for CNMSSM, with an additional one: 

• Ah: by default. 

NMSSMTools allows also for other optional parameters such as //, B', £s an d £f- 

With NMSSMTools 2.2.0, running the program with: 
. /ngmsb. x le5 2el0 5 20 0.1 1 -1000 

returns 

delta0=8. 104e-02 
BR_bsgamma=3 . 089e-04 
BR_Btaunu=l . 094e-04 
Rtaunu=9.928e-01 
BR_Kmunu/BR_pimunu=6 . 454e-01 
R123=1.000e+00 
BR_BDtaunu=6 . 968e-03 
BR_BDtaunu/BR_BDenu=2 . 970e-01 
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BR_Bsmumu=l . 096e-09 
BR_Dstaunu=4 . 818e-02 
BR_Dsmunu=4 . 975e-03 
a_muon=8.081e-10 
collider_excluded=l 
theory_excluded=l 

4.9 NNUHM inputs 

The program nnuhm . x computes the observables using the NNUHM parameters generated 
by NMSSMTools. The necessary arguments to this program are the same as for CNMSSM, 
with two additional ones: 

• tuq: universal scalar mass at GUT scale, 

• mi/2- universal gaugino mass at GUT scale, 

• Aq: trilinear soft breaking parameter at GUT scale, 

• tan/?: ratio of the two Higgs vacuum expectation values, 

• A: cubic Higgs coupling, 

• Mh d ■ down Higgs mass parameter at GUT scale, 

• Mh v : up Higgs mass parameter at GUT scale. 

Optional arguments can also be given: 

• sign(n): sign of Higgsino mass term, positive by default, 

• A K : trilinear soft breaking parameter at GUT scale, by default A K = Aq, 

• mf t ole : top quark pole mass, by default 172.4 GeV, 

• mb(nib): scale independent b-quark mass, by default 4.2 GeV, 

• a s (Mz): strong coupling constant at the Z-boson mass, by default 0.1176. 

In the absence of arguments, a message will be shown. 

With NMSSMTools 2.2.0, running the program with: 
. /nnuhm. x 500 500 50 0.1 500 500 

returns 

SLHA file generated by NMSSMtools 
delta0=9.954e-02 
BR_bsgamma=2 . 232e-04 
BR_Btaunu=6 . 487e-05 
Rtaunu=5.887e-01 
BR_Kmunu/BR_pimunu=6 . 429e-01 
R123=9.980e-01 
BR_BDtaunu=6 . 286e-03 
BR_BDtaunu/BR_BDenu=2 . 679e-01 
BR_Bsmumu=2 . 108e-07 
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BR_Dstaunu=4 . 794e-02 
BR_Dsmunu=4 . 951e-03 
a_muon=2. 121e-09 
charged_LSP=0 
collider_excluded=l 
theory_excluded=0 



4.10 2HDM inputs 

The program thdm . x computes the observables using the 2HDM parameters generated by 
2HDMC. The necessary arguments to this program are: 

• type: Yukawa type (1-4), 

• tan/3: ratio of the two Higgs vacuum expectation values, 

• mA'- CP-odd Higgs mass. 
Optional arguments can also be given: 

• Ai, ■ ■ ■ , A7: Higgs potential parameters, 

• m\ 2 : Higgs potential parameter alternative to uia- 

By specifying the Higgs potential parameters, it is possible to do the calculations in general 
2HDM. If not specified, the optional arguments are set to the default tree level MSSM-like 
values: 

g 2 + g' 2 g 2 - g' 2 g 2 

Ai — A2 — , A3 — , A4 — — — , A5 — A6 — A7 — , 

m i2 = m A cos sin /? . 
In the absence of the necessary arguments, a message will be shown. 

With 2HDMC 1.0.5, running the program with: 

./thdm.x 4 10 300 

returns 

delta0=7.840e-02 
BR_bsgamma=2 . 469e-04 
BR_Btaunu=l . 103e-04 
Rtaunu=1.001e+00 
BR_Kmunu/BR_pimunu=6 . 347e-01 
R123=1.000e+00 
BR_BDtaunu=6 . 979e-03 
BR_BDtaunu/BR_BDenu=2 . 975e-01 
BR_Bsmumu=3 . 273e-09 
BR_Dstaunu=4 . 818e-02 
BR_Dsmunu=4 . 943e-03 
a_muon=-l . 514e-12 

Using the aforementioned main programs as examples, the user is encouraged to write 
his/her own programs in order to, for example, perform scans in a given supersymmetric 
scenario. 
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5 Results 

We illustrate in this section the constraints on the SUSY parameter space that can be 
obtained using observables calculated with Superlso v3.0. For more extended discussions 
about the constraints obtained using Superlso, see for example [13-19]. In Figures Q] and 
[21 two examples of the obtained constraints in the mSUGRA and NUHM scenarios are 
displayed. The different areas in the figures correspond to the following observables: 

• red region: excluded by the isospin asymmetry, 

• blue region: excluded by the inclusive branching ratio of b — > sy, 

• black hatched region: excluded by the collider mass limits, 

• violet region: excluded by the branching ratio of B s — * 

• gray hatched region: favored by the anomalous magnetic moment of the muon, 

• yellow hatched region: charged LSP, disfavored by cosmology, 

• green region: excluded by the branching ratio of B u — > tv t , 

• orange region: excluded by the branching ratio of B — » D°tv t , 

• cyan region: excluded by the branching ratio of K — > /iv^. 

The allowed interval for each observable is given in Appendix iHl 

In Figure HJ the exclusion regions in the mSUGRA parameter plane (?7ti/2 — mo) f° r 
tan (3 = 50, Aq = and \x > are displayed. One can notice that small values of too 
and mi/2 are disfavored by the observables. The unfilled region in the bottom left corner 
corresponds to points with tachyonic particles. 

In Figure [21 the exclusion zones are displayed in the NUHM parameter plane (rriA — tan 0) 
for m = 500 GeV, to 1/2 = 500 GeV, A = and \x = 500 GeV. Most observables tend to 
disfavor the high tan/3 region in this plane. The white top right triangle corresponds to a 
region where tachyonic particles are encountered. 

6 Conclusion 

Superlso v3.0 features many new additions and improvements as compared to the first 
version of the program. It is now able to compute numerous flavor physics observables - 
as well as the muon anomalous magnetic moment - which have already proven to be very 
useful in the exploration of the MSSM and NMSSM parameter spaces. Investigating the 
indirect constraints has many interesting phenomenological impacts, and can provide us 
with important information before the LHC data. They can also be used as guidelines for 
the LHC direct searches, and will be very valuable for the consistency checks. 

In spite of the numerous changes, the spirit of the program is still based on the simplicity of 
use. The code will continue to incorporate other B physics observables. Also, the extension 
of the program to beyond minimal flavor violation is under development. 
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Figure 1: Constraints in mSUGRA (th-i/2 — m o) parameter plane. For the description of 
the various colored zones see the text. The contours are superimposed in the order given 
in the legend. 
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Figure 2: Constraints in NUHM (tua — tan j3) parameter plane. For the description of the 
various colored zones see the text. The contours are superimposed in the order given in the 
legend. 
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Appendix A QCD coupling 



The a s evolution is expressed as [10]: 
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rif denotes the number of active flavors (i.e. 4 for energies between the charm and the 
bottom masses, 5 between the bottom and the top masses, and 6 beyond the top mass). 
We compute the associated A n/ by requiring continuity of a s . In particular in Super I so, 
A5 is calculated so that a s {Mz) matches the value given in the input SLHA file, and then 
A4 and Aq are calculated if needed by imposing the continuity at the bottom and top mass 
scales respectively, which are also input parameters. 

Appendix B Evolution of quark masses 

We use the following three loop formula to compute the pole mass of quarks [10]: 
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For the MS top mass, we use [10] 
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The running mass of the quarks is given by [20] 
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with n/ the number of active flavors, and /3's given in Eq. ([2]). 
The IS" bottom quark mass is given by [21]: 
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In the above equations Ca = 3, T = 1/2, rij ; = 4 and O^ = 4/3. 

Appendix C Wilson coefficients at matching scale 

The effective Hamiltonian describing the b — > transitions has the following generic struc- 
ture: 

n *ff = E v ps v p<> E c ^ °* . ( 17 ) 

* p=u,c 1=1 

where Gf is the Fermi coupling constant, V^- are elements of the CKM matrix, Oj(/z) are 
the relevant operators and are the corresponding Wilson coefficients evaluated at the 

scale jU. 

The Wilson coefficients are given in the standard operator basis [22]: 

01 = (s^T a P L c)(c^T a P L b) , 

2 = (s lfl P L c)(c^P L b) , 

3 = (*7^)E^) ' 

4 = (s llx T a P L b)Y,{q^T a q) , 

: (^i7M 2 7 w ^6)^(^ 1 7 M2 7 /i3 9) , (18) 
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7 = Y^[sa^(m s P L + m b P R )b]F^ , 
Og = ^[^(m^ + m^T^G^ , 
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where Pr,l = (1± 7s)/2- In this basis, the effective Wilson coefficients are defined as [23]: 
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C.l Standard Model contributions 

We express here the SM contributions to the Wilson coefficients following [24,25]. 
The LO coefficients read: 
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and the NNLO coefficients are: 
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243 81 81 27 


r c(2) 




85 5 2 35 L 5 L 2 
~162 ~ 108^ ~ 108 ~ 36 


c *(2) 




= G\{x t w) , 


c *(2) 




= E[{x tw ) , 


c *(2) 




1 2 

= -^G\{x t w) + —El{xtw) , 






3 1 



- 17L 2 , 



(22) 



where 



and [Xw = 0(Mw) 



™ = • < 23 > 

l - Kit) ■ (24) 
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The necessary functions in Eqs. (|20H22|) are: 



At , . -3x 3 + 2x 2 , 22x 3 
Al(x) = — — lnx + 



2(1 - x) 



153x 2 + 159x - 46 
36(1 - x) 3 



(25) 



, -9x 2 + 16x-4, -7x 3 - 21x 2 + 42x + 4 
Eq{x) = — - lnx + 



6(1 - x) 



36(1 - x) 



(26) 



3x 2 



2(1 -x) 4 



lnx + 



5x 3 - 9x 2 + 30x 
12(1 - x) 3 



(27) 



Af , . 32x 4 + 244x 3 - 160x 2 + 16x / 1 

A\ (x) = ; - A Li 2 1 

1V ; 9(1 -x) 4 V x 

-774x 4 - 2826x 3 + 1994x 2 - 130x + 8 , 

H ; te lnx 

81(1 -x) 5 

-94x 4 - 18665x 3 + 20682x 2 - 9113x + 2006 



243(1 - x) 



(28) 



+ 



-12x 4 - 92x 3 + 56x 2 
3(1 - x) 5 



lnx + 



-68x 4 - 202x 3 - 804x 2 + 794x - 152" 
27(1 -x) 4 _ 



In 



in 



2 ' 



, . . 515x 4 - 614x 3 - 81x 2 - 190x + 40 T / 1 
E[{x) = — — Li 2 (l-- 



54(1 - x) 



-1030x 4 + 435x 3 + 1373x 2 + 1950x - 424 



lnx 



+ 



+ - 



108(1 -x) 5 

-29467x 4 + 45604x 3 - 30237x 2 + 66532x - 10960 
1944(1 - x) 4 

133x 4 - 2758x 3 - 2061x 2 + 11522x - 1652 
324(1 -x) 4 

-1125x 3 + 1685x 2 + 380x - 76 
54(1 -x) 5 



lnx 



In 



in 



2 ' 



(29) 
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as well as: 



* . . 4x 4 - 40x 3 - 41x 2 - x T . / 1 

-144x 4 + 3177x 3 + 3661x 2 + 250x - 32 , 
+ 108(1 - xf lUX 

-247x 4 + 11890x 3 + 31779x 2 - 2966x + 1016 



648(1 - x) 4 



+ 



17x 3 + 31x 2 -35x 4 + 170x 3 + 447x 2 + 338x - 56 
mx H 



(1-x) 



18(1 -x) 



In 



. 10x 4 - 100x 3 + 30x 2 + 160x- 40 T / 1 

= 27(13^ L n 1- x- 

30x 3 - 42x 2 - 332x + 68 , -6x 3 - 293x 2 + 161x + 42 
H T^rr, u mx H 



81(1 -x) 4 
90x 2 - 160x + 40 



+ 



27(1 - x) 4 



lnx + 



81(1 -x) 3 

35x 3 + 105x 2 - 210x - 20 
81(1 -x) 3 



In 



m 



2 ' 



T(x) = -(16x + 8)V4x - 1 Cl 2 ^2arcsin^=^ + ^16x + y^ln: 



112 

+32x + — . 



The integral representations for the functions Li 2 and Cl 2 are as follows: 



Li 2 (z) 



f 

Jo 



dt 



ln(l - t) 



Cl 2 (x) = Im [U 2 (e ix )] = - [ d01n|2sin(0/2)| 

Jo 
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The remaining NNLO coefficients take the form: 



n c{2), , ~ C (3), v . 13763 n%f 814 2 \i 2 w 

C 7 M = C 7 i»w = M w ) + — \n w + — In _ 



(35) 



C 8 c(2 W) = cf\^ w = M w ) + 



16607 lni 4 + 397 ^ 2 ^ 



5832 486 M 2 , 



(36) 



C? 3) ( M pf) 



Cy (3) (//w = m t ) 



+ ln 



m? 



-592x 5 - 22x 4 + 12814x 3 - 6376x 2 + 512x . 



27(x - l) 5 



Li 2 1 



(37) 



-26838a; 5 + 25938x 4 + 627367x 3 - 331956x 2 + 16989a; - 460 , 
+ 729(x - 1)6 lnX 

| 34400x 5 + 276644x 4 - 2668324x 3 + 1694437x 2 - 323354x + 53077 



+ ln 2 ^| 



2187(x - l) 5 
-63x 5 + 532x 4 + 2089x 3 - 1118x 2 



lnx 



+ 



9(x - l) 6 

1186x 5 - 2705x 4 - 24791x 3 - 16099x 2 + 19229x - 2740 



162(x - l) 5 



cfW) 



(38) 



+ ln 



H 2 W f -148x 5 + 1052x 4 - 4811x 3 - 3520x 2 - 61x / 1 
~mj [ 18(x - l) 5 2 V x 



-15984x 5 + 152379x 4 - 1358060x 3 - 1201653x 2 - 74190x + 9188 
+ 1944(x - l) 6 

109669x 5 - 1112675x 4 + 6239377x 3 + 8967623x 2 + 768722x - 42796 



+ 



+ ln 2 ^| 

mf 



11664(x - l) 5 

-139x 4 - 2938x 3 - 2683x 2 , 

mx 



+ 



12(x-l) 6 

1295x 5 - 7009x 4 + 29495x 3 + 64513x 2 + 17458x - 2072 



216(x - l) 5 



As regard to the three- loop quantities C^ 3 \fiw = M\y), C^'d^w = M\y), Cy W \nw = m t) 
t(s) 

and C 8 (nw = m t), we have access to their expansions at x — > 1 and x — > oo. 



c(3). 



^(3) , 
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Denoting z = 1/x and w = 1 — z, the coefficients become: 

C^ 3) (fj, w = M w ) ~ 1.525 -0.1165z + 0.01975zlnz + 0.06283z 2 + 0.005349z 2 In z (39) 

+0.01005z 2 In 2 z - 0.04202z 3 + 0.01535z 3 In z - 0.00329z 3 In 2 z 
+0.002372z 4 - 0.0007910z 4 lnz + 0(z 5 ) , 

C<7 {3) (n w = M W ) ~ 1.432 + 0.06709u; + 0.01257w 2 + 0.004710w 3 + 0.002373w 4 (40) 

+0.001406W 5 +0.0009216w 6 +0.00064730w 7 +0.0004779u> 8 +0(w 9 ) , 

C^inw = M w ) ~ -1.870 + O.lOIOz - 0.1218z In z + 0.1045z 2 - 0.03748z 2 lnz (41) 

+0.01151z 2 In 2 z - 0.01023z 3 + 0.004342z 3 In z + 0.0003031z 3 In 2 z 
-0.001537z 4 + 0.0007532z 4 lnz + 0(z 5 ) , 

Cg (3) (//w = M w ) ~ -1.676 - 0.1179™ - 0.02926w 2 - 0.01297w 3 - 0.007296u> 4 (42) 

-0.004672w 5 - 0.003248w 6 - 0.002389u> 7 - 0.001831u> 8 + 0(w 9 ) , 

(fi w = m t ) ~ 12.06 + 12.93z + 3.013zlnz + 96.71z 2 + 52.73z 2 lnz + 147.9z 3 (43) 
+ 187.7z 3 In z - 144.9z 4 + 236.1z 4 In z + 0(z 5 ) , 

cfVw = ra t ) ^ 11-74 + 0.3642u> + 0.1155u> 2 - 0.003145u> 3 - 0.03263u; 4 (44) 
-0.03528u> 5 - 0.03076u> 6 - 0.02504u> 7 - 0.01985u> 8 + 0(w 9 ) , 

Cl i3) (fi w = m t ) ~ -0.8954 - 7.043z- 98.34z 2 -46.21z 2 In z-127.1z 3 (45) 
-181.6z 3 In z + 535.8z 4 - 76.76z 4 In z + 0(z 5 ) , 

Cf i \n w = m t ) ~ -0.6141 - 0.8975^ - 0.03492w 2 + 0.06791w 3 + 0.07966W 4 (46) 
+0.07226u> 5 + 0.06132u> 6 + 0.05096u> 7 + 0.04216u> 8 + 0{w 9 ) . 

For n = 0, 1, 2 and i = 1, ■ ■ ■ ,8, the Wilson coefficients are obtained using: 

^(n) = ^(n) _ ^(n) _ ^ 
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Type 




Add 


All 


I 


cot (3 


cot P 


cot p 


II 


cot P 


— tan/3 


— tan/3 


III 


cot 


— tan/3 


cot /? 


IV 


cot P 


cot P 


— tan p 



Table 1: Yukawa couplings for the four types of 2HDM. U, D and L stand respectively for 
the up-type quarks, the down-type quarks and the leptons. 

C.2 Charged Higgs contributions 

At the Leading Order, the relevant charged Higgs contributions to the Wilson coefficients 
are given by [26]: 

5C?P (fi w ) = (x tH ± ) - A«A 66 if 8 > (x tH ± ) , (48) 



where 



and 



XtH± = ~MiT ' (49) 



m). . s(7 - 5x - 8x 2 ) x 2 3x-2, 

R ' i = -i + -/lnx, 

7 v ' 24(rc-l) 3 4(x-l) 4 

ti), . x(2 + 5x-x 2 ) 3x 2 , . , 

F * W = 8(x-l)» -i(^TF' DI - (50) 

(2) _ g(3 - 5s) x(3x - 2) 

7 lXj " 12(x-l) 2 + 6(x-l)3 mX ' 

(2) _ a(3-s) x 

8 W " 4(x-l) 2 2(x-l) 3 

Xu, A^ are the Yukawa couplings. In super symmetry, they read: 

11 

A« = -j- = • 51 

Aftft tan p 

For the different types of 2HDM, the Yukawa couplings are summarized in Table [TJ 
At the NLO, the charged Higgs contributions can be written in the form [26]: 

<5CfW) = G?(x tH± ) + A?(x tH± )ln-^-^E H (x tH± ), (52) 

m h± y 

<5CfW) = Gf(x^ ± ) + Af(x 4H± )ln^--i^(x^ ± ), (53) 
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4 



4(-3 + 7x-2x») 



3(x- l) 3 



x 



8 - 14a: - 3x 2 , 2 



3(x - iy 



In x 



(54) 



7 - 13x + 2x 2 2(-3 - x + 12x 2 - 2x 3 ) , 
H ; ttt, 1 — — In x 



(x - l) 3 



3(x - l) 4 



;(18 - 37x + 8x 2 ) 
(x^l) 4 " 



Li 2 1 



x 



x(-14 + 23x + 3x 2 ) , 2 
+ — ; — In 2 x 



-50 + 251x - 174x 2 - 192x 3 + 21x 4 , 
ttk In x + 

9(x - l) 5 



(x - l) 5 

797 - 5436x + 7569x 2 - 1202x 3 " 
108(x - l) 4 _ 



21 - 47x + 8x 2 2(-8 + 14x + 3x 2 ) 
(x - 1)3 + (x - l) 4 . 

-31 - 18x + 135x 2 - 14x 3 x(14 - 23x - 3x 2 ) 



(55) 



6(x - l) 4 

-36 + 25x - 17x 2 
2(x- l) 3 

3 - 187x + 12x 2 - 14x 3 



1 

A** A&6 -x 



Lio 1 



+ 



1 



(x - iy 



lnx 



+ 



19 + 17x 2 

m x 



4(x - l) 4 



lnx + 



xj (x - l) 4 

3(143 - 44x + 29x 2 ) 
8(x - l) 3 



(56) 



x(30- 17x + 13x 2 ) T / 1\ x(31 + 17x) l2 
(x-1) 4 Li? ( 1 — \ .J in 2 



-226 + 817x + 1353x 2 + 318x 3 + 42x 4 , 

+ n^ry lnx + 

81 - 16x + 7x 2 19 + 17x 



(x - l) 5 

1130 - 18153x + 7650x 2 - 4451x 3 ' 



216(x - l) 4 



A«Abb-x 



+ x h x 



2(x-l) 3 (x-l) 4lnX 
-38 - 261x + 18x 2 - 7x 3 x(31 + 17s) 



(57) 



+ 



\ 2 



6(x-l) 4 ' (x-1) 5 

x(16 - 29x + 7x 2 ) x(3x - 2) 



lnx 



36(x - l) 3 



+ 



6(x - 1)' 



lnx 



(58) 
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C.3 Supersymmetric contributions 

The relevant chargino contributions are given by [27,28]: 

2 2 ( n 



«?,w = -EE]|ir«i J #^H» f )+ r «n,^(^)} <») 

fc=i i=i L t fc Xi 



2 M W / \ i "p* "p/ M w „(3) 

i V J 2 if, p^pWrr +i+f* f' ^f {3) (t + )\ 

+ Z^i 3 l il *l m 2 ^7,8l a; g a2X ±J+ i li i li-— ^7,8 ^g 12X ±J f » 
i=l I <?12 X, J 

where f/rrij, and ra g - 12 is the common mass of the up and charm squarks, that we 

(3) 

consider identical (m^ 12 ~ ~ mg). The functions iCg(x) take the form 



2 



_(3l. . (5 — 7x) x(3a; — 2) 



F "'M = 2^-(^TF lnl ' (60> 



and Fy 8 (x) are given in Eq. ([50]) . Moreover 

\/2M w sin /3 



lJ >/2cos/?(l + ef tan 0) ' 



(61) 



where U and V are the chargino mixing matrices, following the diagonalizing convention: 

( M 2 M W y/2wn$\ jr _ x 



U (m w ?2^ " (62) 



and D q is the squark q mixing matrix such as: 

D , = (cose, -sin^A 

ysmt^q COSOg J 

€f), which will be given below, is a two loop SUSY correction, with effects enhanced by 
factors of tan (5. 

Tij and V- are obtained from T{j and V- by replacing the matrix D by the unity ma- 
trix. In the following, we adopt the notations: 

COS 9q = D\ x = Cq , sill 9g = D\ x = Sq . (64) 

The leading tan (5 corrections are contained in the following expressions for e' h and e' [29], 
which are evaluated at a typical SUSY scale, fj, s . can be split into two parts: 

e b = e + e 2 , (65) 
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with 



2a s (jK s ) A h / tan/3 - \i 



3tt 



a{Mz)iiM 2 

4 sill 2 0\yTT 



2m? 



H- 



I Ml i? \ 



2 ' 2 
01 01 , 



+ 



2m? 

02 



Ho 



Ml f 



02 02 - 



(66) 



and 



C2 



16 vr 2 



/Lt/ tan (3 — At 



i=i 

a(M z )fiM 2 



H 2 (x hx ±,x i2X± )V i2 



4 sin Ow^ 



c? 



— 2~ 2 — 2~' _ 2" H 2~ ""2 — 2~' — 2~ 

m~ \ mi / m~ \ ra- 

ti \ tl tl / *2 \ *2 *2 / 



(67) 



where A q is the trilinear coupling of the quark q. y q and y q are the ordinary and supersym- 
metric Yukawa couplings of the quark q respectively. The function H 2 is defined as: 



H 2 (x,y) 



xln x 



+ 



2/ In y 



(l-x)(x-y) (l-y)(y-x) 



(68) 



We neglect the neutralino mixing matrices and we assume that the chargino masses are 
given by \x and M 2 . 



e' b {t) 



2a s (ji s ) A b / tan/? - \x 

37T TO; 



C t C b H l( X h~9> X b 2 ~ 9 ) 



+c?s? H 2 (x h - g , x~ bi -) + s\c{ H 2 (x h§ ,x- b2 - g ) + sfsl H 2 (x hg , x~ bi ~ 



,2 2 



2„2 



+ 



16 7T 2 



i=l 



m v o 



(69) 



+c?s? H 2 (x hx o,x~ b2X o) + sf eg H 2 (x ilX o,x biX o) + sfsf H 2 (x ilX o, x^) 
, a(M z )nM 2 



N i3 



4 sin 2 ^vkti" 
o2 / 



^2 



2 

mf 

01 



2 ' 2 

mf mf 

01 Ol , 



s ~6 „ m 
H — ti 2 



2 

m~ 

02 



2 ' 2 
02 02 , 



cj H Ml _sj_ H Ml 

' : 2 \ m? ' m? I 2m? 2 I m? ' m? 

\ tl tl / (2 \ 12 52 , 



2m? 

51 



In the above equation, N represents the neutralino mixing matrix and n x o the number of 
neutralinos, i.e. four in the MSSM and five in the NMSSM. 
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The last e correction reads: 



4 



2a s (/i s ) /i + A t /tan/3 



3vr 



+ 



16 7T 2 



/i/ tan/3 



c 2 l H 2 {x i2 - g ,Xs-g) + s|-f^2(a 
C I4 H ^ x hx^H 2X °) + 4 s l H ^ x hx°> x b ix ^ 



(70) 



+ 4 c f ^(x t - 2X o, ^ o) + sfag H 2 (x hx o,x~ b o 



The SM and charged Higgs contributions at the //jy scale are affected by e^, e' b and eg as 
the following: 



SC. 



(SM, tan 0) 



7,8 



SC. 



(H,tan P) 
7,8 



1 + tan (3 



F 7,i(%tw) , 



{e' + e b ) tan/? (2) 
l + e 6 tan/3 7 ' 8 



(71) 

(72) 



Higher order charged Higgs contributions are expressed as [30]: 

£261 tan 2 /? 



(H,ta,n 2 P) 



^ 8 ^ W> (l + e tan/3)(l + e;) tan/?)" 7 
where F^(x) are given in Eq. ([50]) and e' x reads [31]: 

1 



(73) 



p(2). 



16vr 2 



-"2 I —j — I — 9 tl 2 



2 ' 2 



2 ' 2 



(74) 



Finally, we add the neutral Higgs contributions [32,33]. In the MSSM, they are expressed 



as 



xl*V(0) f„ ^ ° 7(8) 
6 U 7(8)^W) - -^r 



e-2 



72 cos 2 /?(l + eotan/3)(l + e fe tan/3) J 



x £ ^1(4* -4* tan /?) (4 + » 



(75) 



where 07 = 1, as = —3 and, for S = (/t°, H°, A ), 

= (— sin a, cos a, i sin 0) , x u = (cos a, sin a, — i cos /?) . 



(76) 
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In the NMSSM, they are generalized to 



W(8)^W/j - 72 cos 2 /3(1 + eotan/3)(1 + ebtan/?) 2 K") 



3 —2 

m 



-i(Ug* + U%* tan/3)(^ + U«e b ) , 



1 m i 



°7(8)W - 72 cos 2 /S (l + eotari/ 3)(l + e6tari/ 3)2 ^ 



2 —2 

m 



£ -4<P£ + C// 2 * tan /3)(E$ + C#e 6 ) 



x 

where i = F°, J3§ and j = a° v A% 

The CP-even and CP-odd Higgs mixing matrices are respectively [34] 

((cos 9h — sin 6hv6+/x)/ tan (3 cos6h —sinOn \ 

(sin 9 H + cos 9 H v5 + /x) sm6 H cos8 H • (79) 

1 — 1/tan/J — t><5+/xtan/? / 

and 

(cos 0,4 sin /3 cos 6*^ cos /? sin 0^ \ 
-sin6»Asin/3 - sin 6 A cos cos6 A . (80) 
— cos (3 sin /3 J 

where 6h and 9a are the mixing angles, x is the scalar VEV, v ~ 246 GeV and 

\/2^ A + 2kx 

o+ = — 1= 

\/2A\ + kx 

The complete Wilson coefficients Cj are obtained by adding the different contributions. 



(81) 



Appendix D Renormalization group equations 
D.l RGE in the standard operator basis 

The transformations of the Wilson coefficients from the matching scale nw to the scale ^ 
in the standard operator basis (Eq. (|18p ) are given by [35]: 

C^ cS (fi b ) = U^C^{n w ) : (82) 

C (1)cfr (w,) = r)[uMcW eS ( m ) + uWc(°> s (» w )] , (83) 

(7(2)eff (/ife) = ^r^(0)^(2)eff (Mw)+ ^(l)^(l)eff (m) + ^(2)^(0)eff (w) l ; (g4) 
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i 


1 


2 


3 


4 


5 


6 


7 


8 


CLi 


14/23 


16/23 


6/23 


-12/23 


0.4086 


-0.4230 


-0.8994 


0.1456 



Table 2: The numbers a^. 



where r\ = a s (fj,w)/a s (/j,j ) ) and C = {C\, • • • , C$}. The C/"( n ) matrix elements read 

n 8 



u, 



kl 



EE™^V rt ' 



(85) 



j=0 i=l 



The powers aj are given in Table The m^- relevant in our calculations for the are 
given in Tables [3HHJ 

D.2 RGE in the traditional operator basis 

The operators in the traditional basis can be expressed as [23] 



Pi 
Pi 
P 3 

Pi 



<? 

(^7X)E(^2): 
<? 

(«i7^2)E(^S)> 

(5S7^)E(9S^«S)- 



(86) 



The operators Pj and Ps are identical to O7 and Os; an d the C7 and Cs coefficients are 
therefore indistinguishable in both bases. Moreover, the LO Wilson coefficients also coincide 
in both bases. The effective coefficients are nevertheless different [23,36]: 



Cj(/x), for i = 1, ...,6 , 

6 

CV(aO + ^VjCjin), for i = 7 , 



6 



(87) 



Cs(m) + ^2zjCj((i), for i = 8 , 



where y = (0, 0, 0, 0, — |, — 1) and 2 = (0,0,0,0,1,0). In the following, the Wilson coeffi- 
cients at scale /ijy are the non-effective coefficients given in the previous section. 
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VJ . VJt: sJ i 


0092 
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—0 0741 
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0144 
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n 
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VJ 


031 7 
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(00) 
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0.0240 


0.0113 


(00) 
m &2i 


0.8623 











-0.9135 


0.0873 


-0.0571 


0.0209 


(00) 
m 84i 


2.1399 











-2.6788 


0.2318 


0.3741 


-0.0670 


(00) 
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Table 3: Values of the m^' relevant for U^' [35]. 
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-0.1424 
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Table 4: Values of the m^ J relevant for [/^ [35]. 
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i 


1 


2 


3 


4 


5 


6 


7 


8 


(20) 
m 72i 


-212.4136 


167.6577 


5.7465 


-3.7262 


28.8574 


-2.1262 


2.2903 


0.1462 


(21) 
m 72i 


-74.7681 


58.5182 


-13.4731 


3.0791 


7.0744 


2.8757 


3.5962 


-1.2982 


(22) 
m 72i 


31.4443 


-18.1165 


23.2117 


13.2771 


-19.8699 


4.0279 
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2.6149 


(20) 
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(22) 
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15.4051 
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(21) 
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(22) 
m 78i 
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Table 5: Values of the m kli relevant for U kl 



(2) 



[35]. 



The Wilson coefficients at LO relevant for the isospin asymmetry calculation are: 

8 

cfV) = cf\ m )Yk^ (j = l,..., 6), 



/ > '"ji 

i=l 



a 



(0)e// 
7 



(A*6 



i=l 



^ 0)e// (w) 



with 



i=i 

a s (/XH/) 



and where the numbers kji, hi and ^ are given in Table [6] and the a% in Table [2j 
At NLO, the relevant coefficients read, for j = 1, • • • ,6 

cfW) 



i=3 



15 



Ti + C®0*w)fji + 



cfW) 



15 



(90) 



(91) 



(92) 
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1 


2 


3 


4 


5 


6 


7 
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hi 


2.2996 


-1.0880 


-3/7 


-1/14 


-0.6494 


-0.0380 


-0.0185 


-0.0057 


hi 


0.8623 











-0.9135 


0.0873 


-0.0571 


0.0209 


hi 








1/2 


-1/2 














eii 


























fli 








0.8136 


0.7142 














gu 








1.0197 


2.9524 














hi 








1/2 


1/2 














en 


























hi 








0.8136 


-0.7142 














92i 








1.0197 


-2.9524 














hi 








-1/14 


1/6 


0.0510 


-0.1403 


-0.0113 


0.0054 


em 














0.1494 


-0.3726 


0.0738 


-0.0173 


hi 








-0.0766 


-0.1455 


-0.8848 


0.4137 


-0.0114 


0.1722 


93i 








-0.1457 


-0.9841 


0.2303 


1.4672 


0.0971 


-0.0213 


hi 








-1/14 


-1/6 


0.0984 


0.1214 


0.0156 


0.0026 


e-Ai 














0.2885 


0.3224 


-0.1025 


-0.0084 


hi 








-0.2353 


-0.0397 


0.4920 


-0.2758 


0.0019 


-0.1449 


94i 








-0.1457 


0.9841 


0.4447 


-1.2696 


-0.1349 


-0.0104 


hi 














-0.0397 


0.0117 


-0.0025 


0.0304 


em 














-0.1163 


0.0310 


0.0162 


-0.0975 


hi 








0.0397 


0.0926 


0.7342 


-0.1262 


-0.1209 


-0.1085 


9m 














-0.1792 


-0.1221 


0.0213 


-0.1197 


hi 














0.0335 


0.0239 


-0.0462 


-0.0112 


em 














0.0982 


0.0634 


0.3026 


0.0358 


hi 








-0.1191 


-0.2778 


-0.5544 


0.1915 


-0.2744 


0.3568 


96i 














0.1513 


-0.2497 


0.3983 


0.0440 


en 


4661194 
816831 


8516 
2217 








-1.9043 


-0.1008 


0.1216 


0.0183 


fli 


-17.3023 


8.5027 


4.5508 


0.7519 


2.0040 


0.7476 


-0.5385 


0.0914 


37i 


14.8088 


-10.8090 


-0.8740 


0.4218 


-2.9347 


0.3971 


0.1600 


0.0225 


k 


0.5784 


-0.3921 


-0.1429 


0.0476 


-0.1275 


0.0317 


0.0078 


-0.0031 



Table 6: Useful numbers [23]. 
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and 



C, 



(l)e// 
7 



(^b) 



M„(l) e //, , . 8 / 37 39\ „(l)eff , \ 



(93) 



+ ^ ( e 7iQ (^wO»7 + f7iC { 2°\^w) + 57 



(0), 



i=l 



15 



7] 7] 



( 297664 i6 7164416 m 256868 37 6698884 e 

+ I - V 23 orrrnm- ^? 23 + 1 , noo T/23 



V 14283 
37208 



4761 



357075 ' 14283 ' 357075 
r? i_ r? i)4 0) ( W ) + A^ 1) (^), 



where 



(94) 



The numbers eji, fji, gji and appearing in the above equations are given in Table [6j 

Appendix E Calculation of flavor observables 

We provide here the detailed calculation of all the implemented observables for reference. 
The meson masses, lifetimes and form factors, as well as the CKM matrix elements which 
appear in the following are given in Appendix HI Finally, in Appendix [Hi we provide also 
some suggested limits for the observables which can be used in order to constrain the 
implemented models. 



E.l Branching ratio of B — ► X s ^ 

The decay B — > Xg'j proceeds through electromagnetic penguin loops, involving W boson 
in the Standard Model, in addition to charged Higgs boson, chargino, neutralino and gluino 
loops in supersymmetric models. The contribution of neutralino and gluino loops is negli- 
gible in minimal flavor violating scenarios. 
The branching ratio of B — > X s ^f can be written as [37] 



BR(B -> X s7 ) = BR(B -» X c e^) eX p 



v t * s v tb 



V cb 



60 

— [P(E ) + N(E ) + e em ] , (95) 

TTKj 



with 



C 



Vub 



Vcb 



T[B^ X c ev] 



(96) 



Y[B -> X u ev\ 

P{Eq) and N{Eq) denote respectively the perturbative and non perturbative contributions, 
where Eq is a cut on the photon energy. The numerical values of C, BR(i5 — > -^c^^Oexp 

and 

Eq are given in Appendix [U 
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Following [37], we can expand P(E ) as 



P(Eo) = P(°V 6 ) + 



where 



and 



+ 



4n 



4tt 



P[ 2 \n b ) + P^{E Q ^ b ) + P^{E G ^ b ) + O (a 3 s ( Pb )) 



,(2) 



(2)/ 



(O)eff, x r »(0)eff. 



«J=1 



where the K$ can be written in the form: 



K 



(i) 

i7 



Rer«-i 7 r ff ln(^) 2 + 2^( ( 5), for i < 6 



if, 



(!) 182 8 2 (ojeff, f H b \ A i(l),^ 

77 = — ■ i]~ ^77 111 (^^p ) ' ! '77 ! () ' 



A" 



78 



9 27 2 



m 



K, 



(i) _ 



2(1 + %)<4 1) (<5), 



for i,j / 7 . 



The matrix ^°"> efl and the quantities are given by [38]: 



46 



,(1) 



833 1. , . . . i 40 . 
_--[a(z) + *(*)] + —,*, 



._ + 2[a(z) + &(*)]--,*, 



2392 8vr 32 v 56. 

w + ivS + y x '- a(l) + 2, " l) + 8T ! '' 

_76l_J* _16 I 5 
729 9^3 27 6 V ; 3 W V ; 243 

56680 32tt 128_ 896 

-^r + iv! + ^- l6 ' l(l)+32i,(l) + ^r*' r ' 



5710 16tt 64 



729 9^3 27 

44 8 „ 8 

-7T H 27T , 

9 27 9 ' 



1 44 22Q6 

T o(l) + y 6(1) + 12o(z) + 206(z) - — 



where 



and the constant can be written as 



w, 



M X 



m 



IS 



Xb = / dx / dy / dv xyln(v + a;(l — x)(l — v)(l — v + vy)) ~ —0.1684 , 
Jo Jo Jo 



and 
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Y 





32 
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The small-m c expansions of a(z) and b(z) up to 0(z 4 ) read 
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16 

= "77 



3C 3) + lnz + - (Inz 2 + — (lnz) 3 

2 3 SW V2 4 4 V 7 12 vy 



(109) 



+ — I lnz (bo. zY H (lnzr z" 1 + 

4 3 2 4 V ' 12 1 ; ' 



7 vr 2 3 n , 2 

h 21nz - -(lnz) 2 

6 4 4 V 7 



457 5vr 2 



1 



+ I — - — -lnz (lnz) 2 z 4 + m 

l 216 18 72 6 y ' 1 

+ ( - - lnz + -(lnz) 2 ] z 2 + z 3 + -z 4 

1 2 6 2 V ' 9 



4-y +lnz +(lnz) 2 j - 



+ 0(z 5 (lnz) 2 ) 



6(z) 



-3 + ^ - lnz J z - ^-z 3 / 2 + ( 1 + it 2 - 2 In z - ~(ln z) 2 ) 



1376 137 



(110) 



+ ir 2 In z + 2(ln z) 2 z 3 + , 

1 12 9 18 / V 225 30 



2tt z 



+ ^1 In z + 2(ln z) 2 + J z 4 

3 



+ Z7T 



-z + (1 - 2 In z)z 2 + ( -— + ^lnz ) z 3 + z 4 



+ C(z 5 (lnz) 



where £(3) is the Riemann zeta function given in Eq. (115j) . Defining 

5=1- 2E /ml s , 



(111) 



we can write 



54 \ 3 / 12 m c ^m b 



(112) 
(113) 



The remaining NNLO correction to P(Eq) is composed of two parts: 
with 



(114) 



i,i=l,2,7,8 



(115) 



where the contribution for i,j = 3, 4, 5, 6 are neglected, and 
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K. 



(2)A, 

27 



/3 Re 



L( 2 ) 

2 r2 



a(z) + &(*) - 



290 
~8T 



100 
~8T 



(116) 



if 



(2)/3 
17 



1 ^(2)/3o 

6 27 ' 



(117) 



K, 



(2)A) 
77 



54 27 3 v 7 \9 3 / 3 b 



(118) 



A" 



(2)A) 
78 



f 1256 64 2 32 /188 8 2 > o 2 

i(2)/3o 



(119) 



(120) 



9 ' V 27 27 

= 2(1 + 5^)^(5), 
The small-m c expansion of Rer!^\z) up to 0(z 4 ) leads to 

+6840 In z - 1440tt 2 In z - 2520C(3) In z + 120(ln z) 2 + 100(ln zf 

-30(ln z) 4 ) z - (43 - 12 In 2 - 3 In z ) z 3/2 — (11475 - 380tt 2 

v ' ' 243 v ; 1215 v 

+96tt 4 + 7200C(3) - 1110 In* - 15607r 2 lnz + 1440C(3)lnz +990(lnz) 2 



+260(lnz) 3 - 60(lnz) 4 ) z 2 + 



2240vr 2 
243 



.5/2 



2187 



(62471 - 2424tt 2 



-33264C(3) - 19494 In z - 5047r 2 lnz - 5184(ln zf + 2160(ln zf) z 3 



2464 
6075 



7T 2 // 2 + 



15103841 7912 
+ 



546750 3645 



:^ + 



2368 
"8T 



C(3) + 



147038 
6075 



In z 



+ W3 n2lnz + S (lnz)2 " W^ zf ] 1 + c)(:< ' 2(111 : ■ 



where 



Lh = In 



i <? 
mi 



(121) 



(122) 
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The function (f>$^° (5) takes the form 



where F^ 2 ' n ^ is given by [39] 
F( 2 - nf ) = S n{ 5(l - z 



4 



?(<S)L 6 + 4 C \ 
Jo 



z 2 -3 



Li 2 (l -z) + 



i 2 -3 



6(z-l) v ' 6(z-l) 
6z 3 - 25z 2 - z - 18 



ln(l - z) ln(z) 



-^ln 2 (l -z) 



36z 



, X1 x vr 2 -49 + 38z 2 -55z 
ln(l-z)-(l + z)- + 



(123) 



1 


"ln 2 (l-z)~ 




13 


"ln(l - z) 


( tt 2 85 \ 
+ + (-18 + 72j 


1 


2 


1 - z 


+ 


" 36 


1-z 


1 - z_ 



-I + 



(124) 



where the [ln n (l — x)/(l — x)] + are the plus-distributions defined in the standard way, and 



Snf = 49/24 + vr 2 /8 - 2C(3)/3 « 2.474. The remaining 



<(2)A) 



(6) functions are neglected. 



The p^ rem term is more difficult to calculate, as its analytic expression is only known 
in the limit m c 3> nib/ 2. In this limit, we have 

P 2 (2)rcm (£ , W,) - £ <f 0*) Cj° )eS (M) < )rCm (^o, W>) , (125) 

ij=l,2,7,8 

where 



(2)rem 



22 



7>-(2)rem 
-"■27 



^(2)rem 
-"■28 



= K 



z 



l 27 



218 208 
243 ~ "8T Lz? 



+ o\- 

z 



4736 o 1150 r 1617980 20060 w . 
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1664 



729 



729 



19683 



243 



. L c + 
81 v z 



(126) 



(127) 



+ - p^) *£> + |(1 - Lo)K$ + O (i) , (128) 
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7^(2)rem 
J\ 17 



1 (2)rem / 5 3 \ (1) 1237 

6 27 + 1 16 "4 Ld J K78 W 



(1) iz^/ 232 70 2 
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T^(2)rcm 



with 



and 



L c = In 



He 



Ld = Lb — In z = In 



m c {ix c ) 



4f em ( 



rcm ((5) reads 
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where and F^ na ^ are given by [39]: 



= SJ(l -z) + \ 



ln 3 (l - z) 
l-z 



21 

+ tt 



ln 2 (l -z) 



+ 



425 _ tt 2 _ C(3) \ 
96 6 2 J 
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+ 
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2(z-l) 2 Li 3 (z-l) + 



z 3 - 2z 2 + 2z - 3 



2(z-l) 



ln(l - z) 



-140z 4 + 219z 3 - 124z 2 + 28z + 27z 5 + 9z 6 + z 8 - 6z 7 - 6 

I2z(z - l) 3 



Li 2 (z - 1) 



+ 



2z 3 -9z 2 -2z + ll . , -27z 2 + 8z 6 - 9 + 2lz - 3z 3 + 64z 4 - 46z 5 
— ln(l — z) 



4(2-1) 
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Li 2 (l - z) 



+ 



-17z 2 + 4z + 4z 3 + ll T . M , 2z 3 + 13-9z 2 T . 4z - 4z 2 + 1 + z 3 3 
4(^1) Ll3(1 " Z) W^^ U ' {Z) + 6(z-l) ln ^ " *> 

-140z 4 + 219z 3 - 124z 2 + 28z + 27z 5 + 9z 6 + z 8 - 6z 7 - 6 



12z(z - l) 3 
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12 
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-27z 2 + 8z 6 - 9 + 21z - 3z 3 + 64z 4 - 46z 5 



+ 



+ 



4z 
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12z(z - l) 3 
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7T 
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z 72 ' 
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and 



F (2 ' na) = S na 5(l -z) + 



11 



ln 2 (l - 2) 



1 - z 



V _95_ 
+ ' 12 + 144 



ln(l - 2) 
1 - i 



+ 



/ C(3) 905 17vr 2 
V~4 288 + ~ri 



Az - Az 2 + 1 + 2 3 
4(^1) 



Li 3 



1 - z 



2-z 



+ (z- l) 2 Li 3 (z - 1) 



Lis 



2-z 



2Li, 



2-z 



+ 



C(3) 



+ 



-1402 4 + 219z 3 - 124z 2 + 28z + 27z 5 + 9z 6 + z 8 - 6z 7 - 6 



242(2 - iy 
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4(2-1) 

2(3 — 2) . ,., N 
+ 4 ln ( 1 ~ *) 

4z - 4z 2 + 1 + 
12(2 - 1) 



ln(l - z) 



U 2 (z - 1) + 



(1 + z)(2z A - 29z 3 + 732 2 - 57z + 15) 



4z - 4z 2 + 1 + z 3 
Li 2 (l - z) + LiaU) + 



24(2 - l) 3 
(z-3)z 



ln 3 (2 -z) + 



4(2 - 1) " ' 2 

ln 2 (l - 2) + 



Li 3 (l-z) 



42 - 4z 2 + 1 + z 3 , 2 , . 42 - 4z 2 + 1 + z 3 ir 2 



A(z - 1) 



-140z 4 + 219z 3 -124z 2 + 28z + 27z 5 + 9z 6 + z 8 -6z 7 -6 , 
H — — hi(l — z ) 



+ 



24z(z - l) 3 

(1 + z)(2z 4 - 29z 3 + 73z 2 - 572 + 15) 
24(2 - l) 3 



z-l 24 
ln(2 - 2) 



ln(l - z) m(z) 



ln J (l - z) 



(z + 2)(z 3 - 5z 2 + 9z - 35) . 2 „ , 2 5 - 3z 4 - 3z 3 + 34z 2 - 242 + 3 vr 2 
-hr(l-z) 



48 



144 



+ 



[Z Z + 3) 24 + 



vr 2 6z 5 + 72 - 392z 3 + 51z 4 + 219z 2 + 92z 



z 3 - 10z 2 + 6z + 7 
8(^-1) 



C(3) + 



1442(z - 1) 

12z 4 - 754z 3 + 11912 2 + 264z - 761 
^1) 



ln(l - 2) + 



(139) 



where S a = 1.216, S na = —4.795 and Lis(x) = f dyLi2(y)/y. 



«x in Eq. (|13Tj) is denned as 



i 4 V = -n 



The missing Kj, and <j>ij are neglected. 



To estimate the value of p^ rcm for m c < rrib/2, it is necessary to use an interpolation 
method. We consider the linear combination: 



(140) 



P 2 (2)rem (z) = Xl [\r^(z)\ 2 -\r^(0)\ 2 ]+x 2 Re[r^(z)-r^(0)} 



.W/-.M2 |Jl)^nM2 



(2), 



.(2), 



+X3 Re (2) - rf° (0)] + x A z-f- Re (2) + x 5 



(141) 
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i 


1 


2 


3 


4 


5 


6 


7 


8 


di 


1.4107 


-0.8380 


-0.4286 


-0.0714 


-0.6494 


-0.0380 


-0.0185 


-0.0057 



Table 7: Useful numbers for K^ 0) [40]. 



In this equation, x±, ■ ■ ■ , £5 are constants to be determined. Noting that P2 2 ^ rcm (0) = X5, 
we make two different assumptions: 

(a) x 5 = 0, 

(b) x 5 = ~P[ 2 \z = 0) - P 3 (2) (z = 0). 



To determine the four other Xi, we impose that for z 3> 1, Eqs. (|125p and (|14ip coin- 
cide, and by matching the different terms in both equations, the x% can be worked out. In 
particular, we can show that 



oD o 



1 



1 



(142) 



and 



X2 



c (o )cS ( 4019 c (o) cfr _ 1184 c ( 0)efr _ 4C (o )efr + 4 c (o )cfr \ (M3) 
y 486 81 3 J 

Finally, x% and X4 can be determined by choosing two large values for z and requiring 
matching between Eqs. (|125h and (|14ip for both values. 

The two possible determinations of X5 lead to different results, and we therefore compute 
the branching ratio in both cases (a) and (b), and then we give the average value as output, 
as advised in [37]. 



The non-perturbative correction N(Eq) reads [40] 

N(E ) = -1 (#W + rJf<°>) + ^ + ... , 

where r = ffij(/i(y)/mj , A2 ~ (ti^, — m|)/4 is given in Appendix[U and 

i=l 

with dj given in Table El and aj in Table [21 and 



(144) 



(145) 



(0) 



(146) 



In Eqs. (gMD and ([H5|> . 77 = a s (/z w )/a s (/i&). 

The electromagnetic correction e em can be written as [29]: 



a 



Ots{Hb) 



,(0), 



(em). 



(147) 
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where 



and 



with 



kk'(^) = ^(v -1) = — ^— In— , (148) 



/-»(em)/ s /32 _A 40 _JL , 88 w\ (o), . . . 

C$ = I — rj 23 -gg77 23 l^lHq'W (149) 



Mem), , 32 _jl , 32 _x , 640 m 704 is 

Co (iih) = V 23 H 23 H f? 23 ?7 23 , (150) 

8 K ^ b > 575 ' 1449 ' 1449 ' 1725 ' ' K ' 

Mem), , 190 _35 359 _it 4276 _m 350531 _jl 

C 2 (w>) = "377^*? 23 ~~ oTTTF ^ 23 +7^77^ 23 + 



8073 ' 3105 ' 121095 ' 1009125 ' 

2 5956 a 38380 i± 748 ie ,., . 

H 77 23 77 23 -| 7123 7723 . (151 

4347 ' 15525 ' 169533 ' 8625 ' v ; 



Using all the above equations, the inclusive branching ratio of B — > X s j can be obtained. 
E.2 Isospin asymmetry of B — > K*7 

The isospin asymmetry Ao in i? — > if*7 decays arises when the photon is emitted from the 
spectator quark. The contribution to the decay width depends therefore on the charge of 
the spectator quark and is different for charged and neutral B meson decays: 

TjB^K^-TjB^^K*^) 

^ 0± r(£° R*o 7 ) + r(5± k*± 7 ) ' { ' 

which can be written as [41]: 

A = Re(b d - b u ) , (153) 
where the spectator dependent coefficients b q take the form: 

127r 2 f B Q q ( f^, T ^ i f K *m K * 



In the same way as for 6 — > 57 branching ratio, the SUSY contributions induced by charged 
Higgs and chargino loops must be taken into account for the calculation of isospin symmetry 
breaking. 

The functions K\ and K<i q can be written in terms of the Wilson coefficients C% in the 
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traditional basis (see Appendix ID.2|) at scale u, b [41]: 
K X = 



C 2 {Hb) (I lr> — + - ] t\ - (,'iA.r,!,) r 



lln^ + |)F ± -G ± (x c6 ) 
6 fib 6 J 



(155) 



A' 



2q 



^ ' c M + «M) 5 qu + f c 4(M6) + C3( ^ 



iV 



(156) 



+ 



N 4vr 

where iV = 3 and Cf = 4/3 are color factors, and: 



C, {flll )[t]n^ + l-H ± (x db )) +-r 2 
3 /X6 3 



''2 



3 C 3 (/i 6 ) + 3 n/(C4(/x 6 ) + C 6 ( M6 )) - 8(iVC 6 (/i 6 ) + C 5 (m&)) 



F ± ln^ + ... 



/if) 



44 4 

-3-C3GU&) - -n/(C 4 (//b) + C 6 (// 6 )) 



ln^ + .... 

Mo 



(157) 



Here the number of flavors rif = 5, and //o = 0{mb) is an arbitrary normalization scale. 
The coefficient a? reads [42]: 



»'7</*"-) C 7 (//„) + as{lJ A b } CF {C 2 ^ b )G 2 (x cb ) +Cs(fi b )G 8 ) 



+ 



47T 

as{ ^ CF \ C 2 (p h )H 2 {x cb ) + C 8 ( Mfc )H8) 



where fih = yfKhjlb is the spectator scale, and 

n i \ 104 1 to , / \ 
G 2 (x ch ) = — —In h g 2 {x cb ) 

27 nib 



G 8 = - In h 58 , 

3 mj 



with 



11 2tH 2ivr 



+ 



3 ' 



02 0*0 



48 + 30ivr - 5vr 2 - 2ivr 3 - 36((3) + (36 + 6in - 9vr 2 ) In: 



+ (3 + 6i7r) ln 2 x + ln 3 x 



2 

+ 9 



18 + 2vr 2 - 2ivr 3 + (12 - 6vr 2 ) In x + 6ivr ln 2 x + In 3 



+ 



1 

27 



9 + 112zvr - 14vr 2 + (182 - 48ivr) lnx - 1261n 2 x 



833 20ivr 8vr 2 
162 ~ ~W + ~~9~ X 



3/2 



(158) 



(159) 
(160) 



(161) 
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fry 

where £(3) is given in Eq. (|15|) and x c b = — % . The function H2(x) in Eq. (|158|) is defined 

m b 

as: 



^ 2 (x) 



2^ 2 /b/ 



A'" 



3N T B->K* m 2 B J Q - £ 
where /i(u, cc) is the hard-scattering function: 

/ \ / 



dv h(v, x)$±(v) , 



(162) 



o 



h(u,x) = — ^ 



Lio 



V 1 



u — 4x + ie 



Li<; 



u 



7 



1 + 



u — Ax + ie 



u 



7 



2 



(163) 



and Li2 is the usual dilogarithm function given in Eq. (|33|) . 

<£_!_ is the light-cone wave function with transverse polarization, which can be written in 
the form [43]: 



l + 3aiC + ai^(5f - 1) 



(164) 



where u = 1 — u and £ = 2u — 1, and <I>bi is the distribution amplitude of the B meson 
involved in the leading- twist projection. Finally: 



4vr 2 f B fk* f\$Bi(0 f\ *l( 



3.V 7 , /; ' 



(165) 



The first negative moment of <&bi can De parametrized by the quantity A# such as 

*$Bi(0 m B 



(166) 



The convolution integrals of the hard-scattering kernels with the meson distribution ampli- 
tudes are as follows: 



3x ' 



G±(s c ) 



I dx G(s c ,x) , 

Jo 6x 



dxU^x)- 9 ^ 



(167) 



G(s c ,x) 



dx 4>±(x) 



l + x 

3x 2 



with s c = (m c /mb) 2 , and 



G(s, x) = —4 / du uu ln(s — uux — ie) , 



(168) 
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and the Gegenbauer momenta [43]: 



6uu 



4 z 3 

+6 <5 + (3uu + « In u + u In u) + 6 8- (u In -u — u In 

^(i + a + «!^ 3 + (^I + 5c 3 A ) (3e 2 -i) 



(169) 



+ 



, , , , 3 - 



V 112 



+- 5+ (2 + In u + In u) + - 5_ (2£ + In u - In u) . 



(170) 



To compute X±, the parameter X = ln(ms/A^) (l + ^e J</3 ) is introduced to parametrize the 
logarithmically divergent integral dx/(l — x). g < 1 and the phase </? are arbitrary, and 
Aft ~ 0.5 GeV is a typical hadronic scale. The remaining parameters are given in Appendix^ 

Superlso first computes numerically all the integrals and the Wilson coefficients, and then 
calculates the isospin asymmetry of B — > AT* 7 using all the above equations. 

E.3 Branching ratio of B u — > rz/ T 

The purely leptonic decay B u — * tv t occurs via W + and H + mediated annihilation pro- 
cesses. This decay is helicity suppressed in the SM, but there is no such suppression for 
the charged Higgs exchange at high tan /3, and the two contributions can therefore be of 
similar magnitudes. This decay is thus very sensitive to charged Higgs boson and provide 
important constraints. 



The branching ratio of B u — > tv t in supersymmetry is given by [44] 



BR(£ U -» tu t ) 



G F f B \V u bY 
8vr 



T B m B m T 1 £ 



m 



2 \ 2 



m 



B 



m 



B 



m 



tan 2 /? 



n 2 



1 + eo tan (3 



(171) 



where eo is given in Eq. (|66p . and is the B^ meson lifetime which is given in Appendix U 
together with the other constants in this equation. 



The following ratio is usually considered to express the new physics contributions: 



R 



MSSM 



BR{B U tv t ) 



MSSM 



BR(S U -> tv t )sm 
which is also implemented in Superlso. 

In the 2HDM, Eq. (fTTTj) takes the form 
G 2 F f B \V ub \ 2 



m 



B 



m 



H+ 



tan 2 



-1 2 



1 + eo tan (3 



(172) 



BR{B U -» r^ T ) 



8tt 



T B m B m T 1 f- 



2 \ 2 



??17 



1 



???. 



H+ 



XbbK 



(173) 



where the Yukawa couplings Xbb, A Tr can be found in Table Q] for the four types of 2HDM 
Yukawa sectors. 
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E.4 Branching ratio of B — > Dtv t 



The semileptonic decay B — ► Dtv t is similar to 5 U — > riv The SM helicity suppression 
here occurs only near the kinematic endpoint. The branching ratio of B — > Dtv t on the 
other hand is about 50 times larger that the branching ratio of B u — > rz^ T in the SM. 

In supersymmetry, the partial rate of the transition B — * Dive (where f = e,/i or t) 
can be written in function of w as [45] 



dT(B -* Dive) 
dw 



G 2 F \V cb \ 2 m% 

1927T 3 



1 



pv(w) 



l-t(w) 



(174) 



m 



B 



m b 



tan 2 /? 



(m b — m c )m 2 H+ 1 + eg tan f3 



Ps(w) 



where w is a kinematic variable defined as: 

l + (m D /m B ) 2 -(p B -p D ) 2 /m 2 B 



w 



2ni£)/mB 



(175) 



with pd and ps the meson four-momenta, and t(w) = m B + m 2 D — 2w; mn tub- Again, eo is 
given in Eq. (|66|) . and 

t(w) = m B + rnjj — 2wm B) 'm B ■ (176) 



In general 2HDM, Eq. (|174p is replaced by 
dT(B^D£u e ) _ G 2 F \V cb \ 2 m% 



dw 



192vr 3 



mi 



1 - *(«;) 



(177) 



mfeAfcfc - m c \ c 
(m b — m c )m 2 



H+ 



ps(w) 



where the Yukawa couplings X bb , Xee can be found in Table CD for the four types of 2HDM 
Yukawa sectors. 

The vector and scalar Dalitz density contributions read [45] 



pv(w) 



p s (w) 



4 fi+^Vf— Vc- 2 - 1 ) 1 T 1 

V m B J \m B J \ 



my 



t(w) 



1 + 



2JR 1 °^ ' < 178 > 



3 m 



B 



1 + 



in: 



1 + W 2 

AM 



2t{w) V 2t(w) J l-w 
where G(w) and A(w) are hadronic form factors. G(w) can be parametrized as 
G(w) = x [1 - 8p 2 z(w) + (51p 2 - 10)z(w) 2 - (252p 2 - 84)zM 3 ] , 

with 



z(w) 



VwTT-V2 



(179) 

(180) 
(181) 



59 



and A(w) is [46] 

A(u>) = 0.46 ± 0.02 . (182) 

The parameters G(l) and p 2 are given in Appendix [U Integrating Eq. (j!74H over w leads 
as a result to the value of the branching ratio. 



The following ratio 

imv ~ BR(B - D0 eVe ) (183) 

is also considered in order to reduce some of the theoretical uncertainties. It can be calcu- 
lated using Eq. ^JM - 



E.5 Branching ratio of B s — > /i + /i 

The rare decay B s — ► p + fJ-~ proceeds via Z° penguin and box diagrams in the SM, and the 
branching ratio is therefore highly suppressed. In supersymmetry, for large values of tan (3 
this decay can receive large contributions from neutral Higgs bosons in chargino, charged 
Higgs and VF-mediated penguins. 



The relevant operators for this decay are 

„2 



do 
Os 
Op 



P R b a ){U) 



16vr 2 



16vr 2 



m b (s a P R b Q )(lj 5 l) 



where a is a color index and Pl,r = (1 =F 7s)/2- 

The branching fraction for BR(S S — ► fj, + p~) is given by [47,48] 



BR(B S -> fi+fjT 



fLTB s MZ\V tb V t * s \\ll 



64vr 3 



Am 2 , 
v_ 

Ml 



4mJ 
ML 



M 2 Bs \C s \ 2 + 



C P M Bs -2C A 



m. 



(184) 



(185) 



where f Bs is the B s decay constant, M Bg is the B s meson mass and t Bs is the B s mean life, 
all given in Appendix HI 



Ca encloses the SM contributions [48]: 



1.033 / mt(mt) x ' ""' 
sin 2 6 W Vl70 GeV 



CA = ^-li^Sh) , (186) 



which includes both LO and NLO QCD corrections. 
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E.5.1 Supersymmetry 

Cs and Cp include the SUSY loop effects due to Higgs boson contributions, box and penguin 
SUSY diagrams and counter-terms. The charged Higgs contribution is given by [47] 



m^mf tan 2 (5 



(m 2 H ± 



Am H± My^ sin 9w \ mf 



(187) 



The SUSY contributions are expressed as [47] 



C 



Box 
S,P 



m M tan 2 (3 



EE E 



2M W ^ ^ ^ m 2 

i,3=l a=l k,m,n=l 



jaiU^V*! ± ■ z p-Ui 2 V j i )D 1 (x ki ,y ai ,z j i) 



m xf 



(188) 



2 a 2 6 3 

r Pen g _ , nytan p V P a (T^k U- n 



i,j=i a,b=l k,m,n=± 



(189) 



M w [ VajU%V* x ± — l -U i2 Vji \D 2 (yaj,Zij)5 ab 5k m 



[/x*(r [/ «) afc (r c/ ^) fcfe ± M (r^) afe (r l/ «t) fcfe ] jD2 ( 2/a . ;2/fe .) 5 



and the counter term contribution is 

tan 3 f3 



C Count 
S,P 



=F 



EE E h x ±^3(y«)^2(r^) am r^ 



i=l a =l m,n=l 



(190) 



where is the mass of the CP-odd Higgs in the MSSM or the mass of the heavy CP-odd 
Higgs in the NMSSM, m qk = (m u ,m c ,m t ) and 



pa 



1 



2y/2wi 2 w V tb V* s 



VmbKs ^ Mw Va(T u ^) r 



with 



2 

mf. 

2 



JJa 



2 

m u 


2 

x\ 


2 ' 




m + 


m z + 







m qn V t2 (T u ^) na ] ; (191) 



(192) 



Uij and Vij refer to the chargino mixing matrix elements, and Ry to the sneutrino mixing 
matrix. The functions D\{x,y, z), D 2 (x,y) and D^{x) are defined as 



Di(x,y,z) 



xlnx 



(1 — x){x — y)(x — z) 
xlnx 



D 2 (x,y) 



(1 - x)(x - y) 



+ (x <-> y) + (x <-> z) , 
+ (x ^ y) , 



(193) 
(194) 
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D 3 (x) 



xlnx 
1 — x 



(195) 



Note that the function D2(x,y) coincides with H2(x,y) in Eq. (|68j) . V L and Y R are given 
by 



U L \T 



{T u L) 



(T Ur ) 



1 

1 

cos% -sin% 

1 

1 

sin% cos0 t - 



(196) 



(197) 



with the stop mixing angle (—it/ 2 ^ 9^ ^ vr/2). 

In the NMSSM, additional contributions arise from the light CP-odd Higgs [49]: 
1 fv5\ 1 



mi tan f3 



imfa \ x J m| s - m\ 



x < 



+ 



S 2 

yEEE [m^D 3 (y ai )Ua(r u ^ am T-„ 



2 6 3 



Xv 



i=l a=l m,n=l 
2 6 3 



V2sm fw iJ=la=ljfc=1 m x± J 

where Xv is related to the particle masses by 

m 2 H± =m 2 A + mw - X 2 v 2 , (199) 
and vS-/x is related to the CP-odd Higgs mixing angle (see Eq. ([80]) ): 



(198) 



7T 

- + 



2 x tan /? 



To take into account high tan/3 corrections, we multiply Cp and Cs by [48,50] 

1 



(200) 



(l + e fe tan/3) 2 



(201) 



where e& is given in Eq. (|65|) . 



Adding all the contributions to the Wilson coefficients, we can then calculate the branching 
ratio. 
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E.5.2 2HDM 

In pure 2HDM without super symmetry, Eq. (|187p does not hold in general and we have to 
calculate separately different contributions to this equation. For this purpose we extend the 
results of [51] to general Yukawa couplings. The most important contributions (for large 
Yukawa couplings) are given in the following. 



The first contribution to Cs and Cp is from box diagrams involving H + and W + : 

Cs = -Cp = -TTTo ^ 2a ( ^ bb ^ ~^ ss ~ 2 — ~^ tt J *W B+( X H+W> x tw) , (202) 



where x H+w = Mfj+/M w , x tw = m 2 (m t )/M w , and 



log y log x 
x — y \y — 1 x — 1 



(203) 



The second contribution comes from penguin diagrams mediated by neutral Higgs bosons 
with H + and W + in the loop: 



Cl 



m. 



Asm 2 6 w V M 2 



sin 2 a cos 2 a 



+ 



h° 



M H0 



X b b~\ -X ss — 2 — -Xu ) Auu P+(x H + w , xtw) , 

m b mb 



m,, 



4M 2 q sm 2 6 w 



m b 



^bb + n^A ss — 2— -Xu ) A w P + (x H + w ,xtw) ) 



(204) 



where a is the Higgs mixing angle. The third contribution originates from penguin diagrams 
mediated by neutral Higgs bosons, involving H + and G + in the loops: 



C% 



m. 



C C P 



4 sin 2 9 W 
X ( X bb + 
m 



sin 2 a {M 2 H+ - M 2 o ) C os 2 a (M 2 H+ - M 



+ 



M 2 H0 



m s 


nit 
- 2— 


m b 


m b 







(205) 



M 2 



4M 2 sin 2 6 W 



M w 



X bb H -X ss — 2 — -X t t ) XuuP+(x H + w ,x t w) 

m b m b 



These results involve the loop function 



y 



x log x y log y 



x — y \ x — 1 
The last contribution is from self-energy diagrams: 

m.. 



Ct 



4 singly V M lo 



sin 2 a cos 2 a 



+ 



M 2 HQ 



x bb + 



m. 



Xss ) Xf^f^ 



m. 



X H+W + [ X bb -\ X ss Xu 

m b 



m b 

P+(x H+w ,x t w) 



(206) 



(207) 
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m. 



4M 2 Q sin 2 6w 



frig 
rrib 



X H+W + ^bb A ss Xtt 



P+(x H+w ,x t w) 



The Yukawa couplings An are given in Table Q] for the four types of 2HDM Yukawa sectors. 
Adding the four contributions, the charged Higgs contribution to the Wilson coefficients 
can be obtained. 

E.6 Branching ratio of K — > {iv^ 

The leptonic kaon decay K —* fiu^ is also very similar to B u — * tv t , and is mediated via 
W + and H + annihilation processes. The charged Higgs contribution is however reduced as 
H + couples to lighter quarks in this case. 

We consider the following ratio in Superlso in order to reduce the theoretical uncertainties 
from fx [52], which reads in supersymmetry 



BR(K -> jiv^) tk 
BR(vr -> 7V 



V,,. 



Vud 



2 p2 / 1 2/2\2 

% / 1 - mym K 



m 



tan 2 /? 



?7i s / 1 + £q tan /3 



(208) 



where 5 em = 0.0070 ± 0.0035 is a long distance electromagnetic correction factor, the ratio 
fx/fn is given in Appendix [J and eo for the second generation of quarks reads: 



e 



2a s jU 
37rm5 



Ho 



m 



<?L d R 



2 ' 2 



where H^ix^y) is given in Eq. (|68p . 

The additional quantity i?£23 [52] is also implemented in Superlso 



(209) 



Rl23 



V us (K i2 ) Kd(0 + -» 0+) 

X 



m 



K + 



nid \ tan 2 j3 
m s J 1 + eo tan /? 



where ti refers to leptonic decays with i particles in the final state, and + 
sponds to nuclear beta decay. 



0+ 



(210) 



corre- 



In general 2HDM, Eq. (pTUj) reads 



Rl23 



m 



AssA[i[i 



(211) 



where the Yukawa couplings X ss , A w can be found in Table [T] for the four types of 2HDM 
Yukawa sectors. 
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E.7 Branching ratio of D s 



The purely leptonic decays D s — > £i/£ are very similar to K — > and proceed via anni- 
hilation of the heavy meson into W + and H + . The charged Higgs boson contribution can 
only suppress the branching ratio and is therefore slightly disfavored. In supersymmetry 
the branching fraction is given by (where I = e,fj, or r) [53,54]: 

/ o \ 2 



BR(D S 



1 



1 + 



m. 



\m H +) 



m s tan 2 (3 
1 + eo tan (5 



(212) 



where td s and fn s are the L>± meson lifetime and decay constant respectively, which are 
given in Appendix [J together with the other constants in this equation, and 



2a s \i 

3 IT 171?) 



Ho 



m L m l R 



1 ' 2 



with H2(x,y) given in Eq. (f68j) . 

In general 2HDM, Eqs. ([2T2j) becomes 

G 2 



BR(D S -► Ivt 



^\VcsV r D miM Ds T Ds 



(213) 



M Ds 



where the Yukawa couplings A cc , . 
Yukawa sectors. 



1 I — ) ■■'*■•<* Xu ■ (214) 

m H + J (m c + m s ) 

can be found in Table [1] for the four types of 2HDM 



Appendix F Muon anomalous magnetic moment 



The magnetic moment of the muon can be written as 

eh 
2m 



M = (i + ^)— , ( 215 ) 



where eh/2m^ is the Dirac moment and the small higher order correction to the tree level 
is called the anomalous magnetic moment = (g — 2)/2. 

F.l Supersymmetry 

Supersymmetry can contribute to the anomaly through chargino-sneutrino and neutralino- 
smuon loops. The one- loop SUSY contributions to a M are [55] 

n x° 2 ( \ 
H° = S2 £ £ -T^^ml 2 + l"£j 2 )*f + 3-Mrt m n? m ]F?{ Xim ) , 



I07T Z z — ' 12mf dm- 

1=1 rn=l V t lrn 



(216) 
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and 



2 r 2m ± ) 

= I l3r(l c 'l 2 + l<£l 3 )*fW + ^Re[«]if (x fe ) \ , (217) 



where n x o is 4 in the MSSM and 5 in the NMSSM. i, m and k are neutralino, smuon and 
chargino mass eigenstate labels respectively, and 

nf m = V^g'NaX^ + y„N i3 X ml , (218) 
nf m = ^(gN^+g'N^X^-y^NisX*^, (219) 

c£ = 2/^2, (220) 
4 = "9^1, (221) 



where y^ = gm^j y2~mw cos (3 is the muon Yukawa coupling, and the X is the smuon 
mixing matrix. Tl 

xp, = m 2 ± / 77i~ as 



mixing matrix. The functions and depend respectively on Xj m = m 2 o/m 2 - and 



X 



k 



FF(x) = t — ^—j [1 - 6x + 3x 2 + 2x 3 - 6x 2 In si , (222) 
(1 — x) 4 L J 

F ?( x ) = 7T- —T* [l-x 2 + 2xlnxl , (223) 
(1 — x) 6 L J 

2 

FP (x) = - —r \2 + 3x- 6x 2 + x 3 + 6x In xl , (224) 

(1 — x) 4 L J 

= 3 ^ [3-4x + x 2 + 2lnxl . (225) 

2(1 — x) J L J 

The anomalous magnetic moment can also receive at one loop contributions from the Higgs 
bosons, which can be large in the NMSSM [56]: 

„ = Mr(f f , (226) 



i^l^^f tan^ /' -f!* , (227) 



7< 



tfaf = -^tan^/ ^ /' ~\ 2 . (228) 



where G M = g/Ay/lMyy, and U H and t/" 4 are respectively the CP-even and CP-odd Higgs 
mixing matrices given in Eqs. (|79p and (|80p . and m/^ and m ai refer respectively to the 
masses of the three CP-even and the two CP-odd Higgs bosons. 
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In addition to these contributions, we also consider the leading logarithm QED correction 
from two-loop evaluation [57]: 



,SI SY _ . si-SV / i 4 « 1ti M SVSY\ 
V 2 loop — a fi, 1 loop I 1 in 



7T m M / 



(229) 



where -Mstjsy is a typical superpartner mass scale. 



The dominant two loop contributions from the photonic Barr-Zee diagrams with physi- 
cal Higgs bosons are as follows [58] 



,(X7-fO 



a 2 ml 

H 1 



a 2A/f2 ■ 2ft E E Re t A ? A :V] fp,s(m 2 + /Ml) 



hi 



(230) 



- 9 % V (iV c Q 2 ) f y y Re[A^A^] f f (m 2 f I Ml ) , (231) 

8ir 2 M 2 sm 2 e w Jri _ M £ f ?*' 1 



f=t,b,f 



where hi and a; stand respectively for (h°, H°) and ^4° in the MSSM, and (h°, H°, H$) and 
(a°, -Al]) in the NMSSM. N c is the color number and Q the electric charge. 

The couplings of the Higgs to muon, charginos and sfermions in the MSSM are given by: 



\[h°,H°,A°] 



A 



{h ,H°,A } 



X [h°,H°] 



[h°,H°] 



sm a cos a 



cos (3 cos (3 
V2M W 



, tan (3 



(232) 



m 



I Uki Vk2 [ cos a, sin a, - cos 0\ 



+Uk2Vki [ - sin a, cos a, - sin /?] j , 



(233) 



2m t 



m 2 sin/3 



j + /x*[sina, -cosa] + At [cos a, sin a] j {D^fD^ , (234) 
< — //[cos a, sin a] + At, [ — sin a, cos a] ((D^YD^ , (235) 



m~ cos/3 

2m T 
m? cos /3 



j — /x* [cos a, sin a] + A T [ — sin a, cos a] j {DJ^DJ^ , (236) 



where C/ and V are the chargino mixing matrices, and Dj is the sfermion / mixing matrix. 



67 



In the NMSSM, these couplings can be generalized as [56]: 



"X-k 



u i2 

cos (3 

E/^tan/3 

^2/1% 

V2M W 



gm 



'X-k 



\U k2 V k2 U% + g {U kl V k2 U% + U k2 V kl Ug) 
W k2 V k2 u£ - g (U kl V k2 cos/3 + U k2 V kl sin/3) u£ 



tk 



2\/2M, 



w 



gm\ 



h t [A t Uli - \{xU£ +v d Ui )]Re[( J D| 1 )*£>* 2 ] 



+ 



+ 



/2 



h 2 t v u UH- g —(v u Ui-v d Ui) 



D{ 2 



hlv u u!!- 3 -l T ^(v u u t H 1 -v d u!i) 



Dt 



ki 



5m S I 



h b [Ai,Ug - \(xUg + v u Ug)]M(Dl i yDl 2 ] 



+ 



+ 



rJf , 3<? 2 + </ 2 



h 2 b v d U% + 9 —{v u Ull-v d U£) 



D 



k:2 



hlv d ui + - v d u£) 



D 



ki 



Tk 



2 ^ L \hr [A T Ui - A (xUg + v u Ui)] ReKDlTD 



k2\ 



+ 



+ 



hlv d U% + 9 -^-{v u Ull-v d uZ) 



\D 



k2 



h 2 T v d Ui + ^— ^ - v d Ui) 



\Di 



ki ? , 



where v u , v d and x are the VEV of H u , H d and S such as 



«2 = 



sin 2 /3 



tan (3 = — . 



The loop integral function fps is given by: 

f 1 1 z(l-z) 2x 

JPS{x) = x / dz— ; In 

Jo z(l - z 



z) — x x y 

with y = \/l — Ax. The other loop functions are related to fps as 

f s (x) = (2x - l)f PS (x) - 2x{2 + logx) , 



f f ~(x) 



X 

2 L 



2 + logx- f PS (x) 



(246) 
(247) 



The contribution from the bosonic electroweak two loop diagrams can be written as [56,60]: 



6a 



bos 



5 Gf mfi ot 
24^2 vr 3 



m. 



CLln W 



+ c 



w 



where 



Ci = 30 



2 \2' 



98 + 9c h L + 23 (1 - 4s 2 w ) 



Cr in the MSSM is given by 



2 r 



^ _ cos 2(3 M% 
L cos/3 

and in the NMSSM is extended to 



cos a cos(a + 0) sin a sin(a + (3) 



+ 



m 



ii 



c\ = cos 2/3 M| 



i=l 



mi 



from which the SM bosonic electroweak two loop contributions have to be deduced: 



5a 



SM 



with 



5 Gf rn^ a 
24^/2 vr 3 

1 



cf M ln 



m 



r SM _ 
Cl ~ 30 L 



2 \ 2 



107 + 23(l-4s^) 



(248) 
(249) 

(250) 
(251) 

(252) 
(253) 



and co and Cq M are neglected. 



F.2 2HDM 

In the 2HDM the dominant contribution can be obtained by generalizing the results of [59]: 



Sat 



2Tr 



f 



(254) 



rn 



H 



k sin(/3 — a) + // cos(/3 — a) k m sin(/3 — a) + cos(/3 — a) f (xfh ) 



k cos(/3 — a) — pf sin(/3 — a) cos(/3 — a) — /Z 1 sin(/3 — a) f (xfn ) 
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where k/ = \/2mf/v, pf = XffK?, Xfx = m 2 /m 2 x , and is the third component of the 
weak isospin, —1/2 for down-type fermions and 1/2 for up-type fermions. Finally, the / 
and g functions are given by 



and 

y{1 ^± . (256) 

x 



s{x) = Sl Ay vil-y)-^ - 



Appendix G LHA file format for 2HDM 

Superlso needs a Les Houches Accord (LHA) inspired input file for the calculations in 
2HDM. This file, in addition to the usual MODSEL, SMINPUTS, GAUGE, MASS and 
ALPHA blocks of SLHA format, needs three additional blocks to specify the Yukawa cou- 
pling matrices for up-type and down-type quarks and for leptons. Also, the 2HDM model 
must be specified in the MODSEL block with the entry followed by a positive integer. 
Such a file can be generated by 2HDMC. An example is given in the following: 



Block MODSEL # Select Model 
10 # 10 = THDM 

Block SMINPUTS # Standard Model inputs 



1 1.27910000e+02 # l/alpha_em(MZ) SM MSbar 

2 1. 16637000e-05 # G Fermi 

3 1. 17600000e-01 # alpha_s(MZ) SM MSbar 

4 9. 11876000e+01 # MZ 

5 4.24680224e+00 # mb(mb) 

6 1.71200000e+02 # mt (pole) 

7 1.77700000e+00 # mtau(pole) 
Block GAUGE # SM Gauge couplings 

1 3.58051564e-01 # g> 

2 6.48408288e-01 # g 

3 1.47780518e+00 # g_3 

Block MASS # Mass spectrum (kinematic masses) 

# PDG Mass 

25 9. 11123204e+01 # Mhl , lightest CP-even Higgs 

35 5.00013710e+02 # Mh2, heaviest CP-even Higgs 

36 4.99999960e+02 # Mh3, CP-odd Higgs 

37 5.06332023e+02 # Mhc 

Block ALPHA # Effective Higgs mixing parameter 
3. 12022092e+00 # alpha 



70 



Block UCOUPL 



1 


1 


0, 


,00000000e+00 


# 


LU_ 


.{1,1} 


1 


2 


0, 


.00000000e+00 


# 


LU_ 


.{1,2} 


1 


3 


0, 


.00000000e+00 


# 


LU_ 


.{1,3} 


2 


1 


0, 


.00000000e+00 


# 


LU_ 


.{2,1} 


2 


2 


2. 


.00000000e-02 


# 


LU_ 


.{2,2} 


2 


3 


0. 


.00000000e+00 


# 


LU_ 


.{2,3} 


3 


1 


0. 


,00000000e+00 


# 


LU_ 


.{3,1} 


3 


2 


0, 


,00000000e+00 


# 


LU_ 


.{3,2} 


3 


3 


2. 


,00000000e-02 


# 


LU_ 


.{3,3} 



Block DCOUPL 



1 


1 


0, 


,00000000e+00 


# 


LD. 


.{1,1} 


1 


2 


0, 


.00000000e+00 


# 


LD. 


.{1,2} 


1 


3 


0, 


.00000000e+00 


# 


LD. 


.{1,3} 


2 


1 


0, 


.00000000e+00 


# 


LD. 


.{2,1} 


2 


2 


-5. 


,00000000e+01 


# 


LD. 


.{2,2} 


2 


3 


0, 


.00000000e+00 


# 


LD. 


.{2,3} 


3 


1 


0, 


.00000000e+00 


# 


LD. 


.{3,1} 


3 


2 


0, 


,00000000e+00 


# 


LD. 


.{3,2} 


3 


3 


-5. 


,00000000e+01 


# 


LD. 


.{3,3} 



Block LCOUPL 



1 


1 


-5. 


,00000000e+01 


# 


LL. 


.{1,1} 


1 


2 


0. 


.00000000e+00 


# 


LL_ 


.{1,2} 


1 


3 


0, 


,00000000e+00 


# 


LL. 


.{1,3} 


2 


1 


0, 


.00000000e+00 


# 


LL. 


.{2,1} 


2 


2 


-5. 


.00000000e+01 


# 


LL. 


.{2,2} 


2 


3 


0, 


,00000000e+00 


# 


LL. 


.{2,3} 


3 


1 


0, 


,00000000e+00 


# 


LL. 


.{3,1} 


3 


2 


0, 


,00000000e+00 


# 


LL. 


.{3,2} 


3 


3 


-5. 


,00000000e+01 


# 


LL. 


.{3,3} 



The UCOUPL, DCOUPL and LCOUPL blocks contain respectively the Yukawa couplings 
^UUi ^dd and If some entries in these blocks are missing, they will be set to 0. The 
Yukawa couplings for the types I-IV are specified in Table [TJ 



71 



Appendix H Suggested limits 



The limits on the masses of Higgs and SUSY particles from direct searches at collid- 
ers are given in Table [5J Some of the limits are subject to auxiliary conditions (see 
[10]) which are also taken into account in the program. These values are encoded in 
src/excluded_masses . c and can be updated by the user if necessary. 

Finally, in Table EJ we present our suggested limits for each observable, which can be 
used to constrain SUSY parameters. 

We would like to stress however that some of the inputs in this table suffer from large un- 
certainties from the determination of CKM matrix elements and/or hadronic parameters. 
The constraints obtained using these observables should therefore not be over-interpreted 
(see [16] for more details). 



Lower bounds on the particle masses in GeV 


h° 


111 


en 


73 


H+ 


79.3 




94 


A 


93.4 


h 


81.9 


X? 


46 


V 


94 


X% 


62.4 


h 


95.7 


xl 


99.9 


h 


89 


xl 


116 


q 


379 


xt 


94 


9 


308 



Table 8: Limit on the masses of Higgs and MSSM particles from direct searches [10]. The 
limit on the h mass includes a 3 GeV intrinsic uncertainty. 



Appendix I Useful parameters 

The masses of quarks and mesons, as well as some other useful parameters such as lifetimes, 
CKM matrix elements or decay constants, are given in Table PTOl 
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Observable 


Combined experimental value 


95% C.L. Bound 


BR(B X sl ) 


(3.52 ± 0.23 ± 0.09) x 10~ 4 [70] 


2.15 x 10~ 4 < BR(6 -► 87) < 4.89 x 10~ 4 




^o.i it z.o) x iu ^aj 


1 7 v 1 PI -2 A ^ on w i n — 2 
— L.I X 1U <. ZAq <. o.y X 1U 


DL\\O u — > TU T ) 
Rtv t 


1.28 ±0.38 (b) 


n QQ v If! - 4 ^ Dp / f> , _,. \ 9 An v in— 4 

0.52 < R TUt < 2.04 


BR(B -► D°ri/ r ) 


(8.6 ±2.4 ±1.1 ±0.6) x 10~ 3 [73] 
0.416 ±0.117 ±0.052 [73] 


2.9 x 10-3 < BR(£ -»• D°tv t ) < 14.2 x 10~ 3 
0.151 < i mv < 0.681 


BR(B S -► M + AO 


< 5.8 x 10- 8 [74] 


BR(B S -► n+fj,-) < 6.6 x 10~ 8 


BR( J FC -► /«/) 

BR(7T — ► /Xl/J 


0.6358 ±0.0011 (c) 
1.004 ± 0.007 [52] 


0.6257 < ^ < 0.6459 

BR(7T — ► [IV) 

0.990 < i? m < 1-018 (d) 


BR(D S -► r^ T ) 
BR(L> S /^) 


(5.7 ±0.4) x 10~ 2 [54] 
5.8 ± 0.4 x 10~ 3 [54] 


4.8 x 10~ 2 < BR(D S -► r^ T ) < 6.6 x 10~ 2 

4.9 x 10~ 3 < BR(D S -► < 6.7 x 10~ 3 


5% 


(2.95 ± 0.88) x 10~ 9 [75] 


1.15 x 10~ 9 < <5a^ < 4.75 x 10~ 9 



Table 9: Suggested limits for the observables implemented in Superlso v3.0. The 95% 
C.L. bounds presented in this table include both the experimental and the theoretical 
uncertainties. 

(a) Value obtained combining the most recent measurement of Babar [71] with the results 
of [10,72]. 

(b) Value deduced from [70]. 

(c) Value obtained combining the results of [10,52]. 

(d) See [16] for a discussion on the uncertainties. 
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Meson masses in MeV [10] 




m K 


m K * 


m D 




m B 




139.6 


493.7 


891.7 


1864.8 


1968.5 


5279.1 


5366.3 





Meson lifetimes in s 


[10] 






TK 


TB 






2.6033 x 10~ 8 


1.2380 x 10~ 8 


1.638 x 10~ i2 


1.425 x 10~ i2 


5.00 x 10~ i3 




Quark masses in GeV and a s [10] 


m b (m b ) 


m c (m c ) 


m s 


[61] 


a s {M z ) 


4 20+ - 17 
4 - zu -0.07 


i 97 +0.07 


n 1 n4 +0 - 026 

U.1U4„ 034 


172.4 ± 1.2 


0.1176 ±0.002 



Hadronic parameters in MeV (fi = 1 GeV) 


[62] 


/a-. [63] 




fB [64] 


As [65] 


fB. [64] 


f Ds [66] 


1.189 ±0.007 


220 ±5 


185 ± 10 


200 ± 20 


460 ± 110 


245 ± 25 


241 ±3 



Meson mass and coupling parameters (fi = 1 GeV) [43] 


S3 


<>3 


<>3 


w l,0 


h 


6- 


0.032 


0.013 


0.024 


-2.1 


0.16 


-0.16 



Meson related parameters (fi = 1 GeV) 


[67] 


[63] 


a\(K*) [67] 


4(K*) [67] 


0.04 ± 0.03 


0.15 ±0.10 


0.03 ± 0.02 


0.11 ±0.09 


G(l) [46] 


P 2 [45] 


A(w) [46] 




1.026 ±0.017 


1.17 ±0.18 


0.46 ± 0.02 


0.31 ±0.04 



CKM parameters 


|Ks| [52] 


\V ub \ [10] 


\Vcb | [10] 


[io] 


0.2253 ± 0.0009 


(3.95 ± 0.35) x 10~ 3 


(4.12 ±0.11) x 10~ 2 


0.97334 ± 0.00023 


IWK<i| [52] 


\V ts V tb /V cb \ 2 [37] 


Re(K* s WK* s K 6 ) t 41 ] 


0.2321 ±0.0015 


0.9676 ± 0.0033 


0.011 ±0.005 



b — > S7 related parameters 


C [68] 


E (GeV) [37] 


A 2 (GeV 2 ) [10,37] 


BR(B X c e^) CX p [69] 


0.58 ±0.01 


1.6 


0.12 


0.1061 ±0.0017 



Table 10: Input parameters. 
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are acknowledged for clarity. We refer to [76] for a fairly complete review of available tools 
in SUSY. 
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